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Mr.  W.  Gordon  Brady  of  CAL  conducted  the  study,  and  Mr.  J.  McHuj 
and,  subsequently,  Mr.  P.  Cancro  administered  the  project  for 
USATRECOM. 

The  authors  extend  their  appreciation  to  Messrs.  C.  Ryan,  J.  Balcer 
and  J.  Nemeth  of  CAL  for  their  significant  contributions  to  the  investi 
gation. 

The  work  reported  herein  is  a  continuation  of  the  initial  effort  under 
Contract  DA  44-177'TC-782  that  was  reported  in  THECOM  Technical 
Report  63*  11,  Theoretical  and  Experimental  Studies  of  Impir.a*ing 
Uniform  Jets,  April  1963. 


iii 


TR  64-42 


TABLE  OF  CONTENTS 


FOREWORD 

LIST  OF  ILLUSTRATIONS 
LIST  OF  SYMBOLS 
SUMMARY 
CONCLUSIONS 
RECOI^ENDATIONS 
TECHNICAL  DISCUSSION 

I  INTRODUCTION 

II  EXPERIMENTAL  PROGRAM 

III  THEORETICAL  ANALYSES 

IV  COMPARISON  OF  THEORY  AND  EXPERIMENT 
AND  DISCUSSION  OF  RESULTS 

BIBLIOGRAPHY 

APPENDIX  I  GENERAL  COMMENTS  ON  THE  FLOW 

VISUALIZATION  PICTURES 

APPENDIX  II  DETAILS  OF  UNIFORM  INVISCID 

IMPINGING  JET  ANALYSIS 

DISTRIBUTION 


IX 

1 

2 

4 

5 

5 

6 
16 

25 

34 

35 

38 

92 


42 


V 


TR  64-42 


1ST  OF  ILLUSTRATIONS 


Figure  Page 


1  Shear  Screen  Assembly  50 

2  Flow  Visualization  Apparatus  51 

3  Flow  in  Plane  of  Symmetry  of  Normally 

Impinging  Uniform  Jet,  H/D  =  2  52 

4  Flow  Visualization  Pictures  of  Normally 
Impinging  Nonuniform  Jet,  H/D  1.  99 

(a)  Flow  in  Plane  of  Symmetry  53 

(b)  Flow  Along  Ground  54 

5  Flow  Visualization  Pictures  of  Tilted 
Impinging  Nonuniform  Jet,  H/D  =  1.99 

(a)  Flow  in  Plane  of  Symmetry  55 

(b)  Flow  Along  Ground  56 

6  Notation  for  Flow  in  Impinging  Jet  57 

7  Velocity  Profiles  in  Free  Nonuniform  Jet 

at  Nozzle  Exit  58 

8  Mean  Velocity  Profiles  in  Free  Nonuniform 

Jet  at  Various  Distances  From  Nozzle  59 

9  Variation  of  Max' mum  Velocity  and  Center* 
line  Velocity  in  Nonuniform  Free  Jet  with 

Distance  from  Nozzle  60 

10  Variation  of  Mass  Flow  in  Uniform  and  Non- 

uniform  Free  Jets  with  Distance  from  Nozzle  61 

11  Mean  Velocity  Profiles  in  Nozzle  Exit  Plane 

for  Various  H/D,  Nonuniform  Jet  62 

12  Mean  Static  Pressure  Profiles  in  Nozzle  Exit 

Plane  for  Various  H/D,  Nonuniform  Jet  63 

13  Ground  Static  Pressure  Distributions,  Normally 
Impinging  Nonuniform  Jet 

(a)  H/D  =  0.49  64 

(b)  H/D  =  0.99  65 

(c)  H/D  =  1.99  66 

(d)  H/D  =  4.  07  67 

V  i 


TR  64-42 


Figure  Page 


14  Ground  Static  Pressure  Distribution,  Non- 

uniform  Jet,  H/D  =  1.  99,  Tilt  Angle  ts  5  degrees  68 

15  Static  Pressure  Distributions  on  Jet  Center- 

line,  Nonuniform  Jet  69 

16  Absolute  Velocity  Distributions  in  Turning 

Region  of  Impinging  Nonuniform  Jet  at  Several 
Constant  Distances  from  Ground,  H/D  =  1.  99  70 

17  Velocity  Profiles  Near  Ground,  Normally 
Impinging  Nonuniform  Jet 

(a)  H/D  =  0.49  71 

(b)  H/D  =  1.99  72 

(c)  H/D  =  4.  07  73 

18  Velocity  Profiles  Near  the  Ground  for  r/D>  1.  0, 

Nonuniform  Jet,  H/D  =  1.  99  74 

19  Boundary- Layer  Profiles  for  r/D>1.0,  Nonuni¬ 
form  Jet,  H/D  =  1.99  75 

?0  Computer  Jet  Boundaries  for  Uniform  Inviscid 

Jet  76 

21  Jet -Centerline  Static  Pressure  Distributions, 

Uniform  Jet  77 

22  Ground-Plane  Static  Pressure  Distributions, 

Uniform  Jet  78 

23  Velocity  Distributions  at  Nossle  Exit,  Uniform  Jet  79 

24  Notation  Used  in  Calculation  of  Radial  Momentum 

Flux  of  Flow  Under  Nonuniform  Impinging  Jet  80 

25  Nondimensional  Radial  Mass  Flow  vs  r/R  for 

Both  Uniform  and  Nonuniform  Jets  81 

26  Variation  of  Flow  Properties  Near  Ground  with 
Radial  Distance  from  Stagnation  Point,  Normally 
Impinging  Uniform  Jet 

(a)  Maximum  Velocity  Near  Ground  82 

(b)  Flow  Thickness  83 


vii 


TR  64-42 


Figure 


Pa^ 


27 

Variation  of  Flow  Properties  Near  Ground 
with  Radial  Distance  from  Stagnation  Point, 
Normally  Impinging  Nonuniform  Jet 

(a)  Maximum  Velocity  Near  Ground 

(b)  Flow  Thickness 

84 

85 

28 

Computed  Absolute  Velocity  Profile  Per¬ 
pendicular  to  Ground  Plane  for  r/D  =  !•  7, 

H/D  =1.0,  Inviscid  Uniform  Jet  Analysis 

86 

29 

Variation  of  U/(J^  with  H/D  lor  a  Uniform  Jet 

87 

30 

Comparison  of  Calculated  Maximum  Velocity 

Near  Ground  at  Start  of  Wall  Jet  under  Uniform 
and  Nominiform  Jets  of  Equal  Thrust  Loading 

88 

31 

Representation  of  Vorticity  Layer  on  Edge  of 
Inviscid  Jet 

89 

32 

Flow  Model  used  in  Inviscid  Analysis  of  Circular 
Impinging  Uniform  Jet 

90 

33 

Breakdown  of  Regions  of  Integration 

91 

li 


viii 


TR  64-42 


USX  OF  SYMB01.S 


D 

H 

Pa 

Pi 

Pi 

R 

U 

Urn 

IV 

Hr 

r 

t 

u 

V 


Diameter  of  Jet  Nozzle  at  Exit  =  1  foot 
Distance  of  Jet  Nozzle  Exit  Frcm  Ground,  feet 
Radial  Momentum  Flux.  Slugs >ft. /(tec. ) 

Atmospheric  Pressure,  psf. 

Static  Pressure,  psf. 


Total  Pressure,  psf. 


Radial  Mass  Flow,  Slugs/ sec. 

2  r 

Avei  age  Velocity  in  Nozzle  ~  •  ^P** 


Mean  Radial  Momentum  Velocity  * 


Reference  Velocity  • 


.  fps. 


—  ,  fps.  (Uniform  Jet  Only) 


Jet  Reference  Velocity  = 


Absolute  Value  of  Velocity,  fps. 

Reference  Dynamic  Pressure  =  Wind  Tunnel  Settling 
Chamber  Pressure,  psf. 


Radial  Distance  Along  Ground  Board  from  Stagnation 
Point,  feet 


Radial  Distance  from  Jet  Axis  in  Nozzle  Exit  Plane,  feet 

Distance  from  Nozzle  Exit,  Positive  Toward  Ground 
Plane,  feet 


Local  Velocity  Parallel  to  Ground  Plane,  fps. 

Maximum  Value  of  U  Attained  in  the  Ground  Flow  at 
Any  One  Radial  Station,  fps.  (See  Figure  6) 

Local  Velocity  Perpendicular  to  Ground  Plane,  fps. 


ix 


TR  64-42 


y 


Velocity  Based  on  Local  Total  Pressure  = 

Distance  Above  G  ound  Board,  feet 

y  When  u.- (see  Figure  6) 

y  When  u  =  ~U^  in  Outer  Ground  Flow  Region 
(see  Figure  6) 


2(Pt -Pa) 


fps. 


2 

cc 

y 


Distance  Above  Ground  Board  Along  Jet  Centerline,  feet 

Form  Parameter  Associated  With  Ground  Flow 
Reynolds  Number  for  Radial  Wall-Jet  (see  Reference  4) 

Angular  Coordinate  About  Jet  Centerline  (  S  is  measured 
clockwise  looking  downstream.  0*0  ta  vertically  upward.  ) 

3 

Axr  Density,  Slugs/ft 

Kinematic  Viscosity,  lx  !  sec 


X 


TR  64-42 


SUMMARY 


The  results  of  an  experimental  investigation  of  the  flow  under  a 
normally  impinging  nonuniform  jet  are  presented.  The  jet  velocity 
profile  was  designed  to  be  representative  of  rotors  and  ducted  fans. 
Primary  emphasis  was  placed  on  determining  the  properties  of  the 
flo'/  near  the  ground.  Some  measurements  of  the  flow  in  the  turning 
region  of  the  jet  are  also  presented  along  with  flow  visualization 
pictures.  The  jet  was  tested  at  distances  from  the  ground  of  4,  Z, 
and  ^  nozzle  diameters.  The  data  are  compared  to  corresponding, 
previously  published  data  obtained  with  a  uniform  jet.  An  approximate 
analysis  which  uses  an  empirical  relation  for  radial  mast  flow  near 
the  ground  is  used  to  calculate  the  properties  of  the  flow  along  the 
ground  at  radii  large  enough  so  that  the  pressure  gradient  is  approxi¬ 
mately  zero.  The  results  of  the  approximate  analysis  are  compared 
to  the  experimental  data  for  both  the  uniform  and  the  nonuniform 
impinging  jets. 

A  method  of  calculating  the  properties  of  the  flow  in  an  inviscid, 
normally  impinging,  uniform  jet  has  been  formulated.  The  formula¬ 
tion  is  applicable  for  all  distances  between  the  jet  nozzle  and  the  ground. 
Solutions  have  been  obtained  for  jets  at  nozzle-to-ground  distances  of 
^  and  1  jet  diameters.  The  mathematical  model  used  was  based  on 
a  vortex-sheet  representation,  and  solutions  were  obtained  by  means 
of  an  iterative  technique  using  an  IBM  704  digital  computer.  Good 
agreement  was  obtained  with  experimental  ground-plane  and  jet- 
centerline  pressure  distributions  .and  with  nozzle-exit  velocity  profiles. 
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CONCLUSIONS 


On  the  basis  of  the  results  presented  both  in  this  report  and  in  Refer¬ 
ence  1,  the  following  conclusions  can  be  drawn  regarding  the  flow  under 
impinging  uniform  and  nonuniform  jets  of  circular  cross  section.  The 
nonuniform  jet  considered  here  is  one  with  a  triangular  velocity  dis¬ 
tribution  similar  to  that  generated  by  rotors  and  ducted  fans. 

1.  The  flow  in  the  turning  region  of  a  normally  impinging  non- 
uniform  jet  with  a  triangular  velocity  distribution  is  sub¬ 
stantially  different  from  that  under  a  uniform  jet.  Under 
the  nouuniform  jet  there  is  a  region  of  upwash  near  the 
center  of  impingement  which  is  not  present  with  a  uniform 
jet. 

?,,  Once  the  jet  has  been  completely  turned  by  the  ground,  the 
flows  under  both  types  of  jets  rapidly  approach  the  form  of 
a  turbulent  radial  wall  jet.  For  the  same  thrust  per  unit 
area^the  maximum  velocity  and  thickness  of  the  flow  in  the 
wall  jet  region  are  almost  identical  for  both  jets  for  values 
of  H/Disl.  At  larger  values  of  H/D,  the  flow  under  the 
nonuniform  jet  is  thicker  and  reaches  a  lower  maximum 
velocity  than  that  under  a  uniform  jet  of  equal  radius  and 
thrust.  These  differences  are  relatively  small. 

3.  The  flow  under  a  nonuniform  jet  with  a  triangular  velocity 

distribution  is  extremely  sensitive  to  small  deviation  angles 
from  the  normal  impingement  condition.  (No  check  of  the 
effect  of  tilt  angle  was  made  during  the  earlier  uniform  jet 
experiments;  but  based  on  the  general  experience  gained 
in  setting  up  conditions  for  normal  flow  impingement  for 
both  types  of  jets,  a  uniform  jet  may  not  exhibit  such 
sensit.vity.  ) 
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4.  An  approximate  analysis  developed  to  predict  the  properties 

of  the  flow  along  the  ground  under  both  uniform  and  nonuniform 
jets  provided  satisfactory  agreement  with  the  experimental 
data  for  sufficiently  large  radii.  (This  analysis)  based  on 
Glauert's  wall-jet  theory  and  an  empirical  curve  for  radial 
mass  flow,  is  strictly  valid  only  in  regions  where  the  static 
pressure  gradient  along  the  ground  is  zero.  ) 

5.  Solutions  for  the  properties  of  the  inviscid  flow  under  uniform 
impinging  jets  with  H/O  =  1.0  and  H/D  =  0.  25  have  been  obtained 
on  a  digital  computer  using  a  vortex-sheet  representation 

for  the  free-jet  boundary.  The  computed  results  for  the 
flow  in  the  turning  region  of  the  impinging  uniform  jet  can 
be  applied  to  the  real  viscous  jet  problem  with  good  accuracy 
in  the  range  H/D^l.  At  larger  values  of  H/D,  viscous  mixing 
significantly  alters  the  real  flow  from  the  inviscid  flow. 
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RECOMMENDATIONS 


On  the  basis  of  the  results  presented  herein  and  in  Reference  1,  the 
following  recomnneudations  are  made. 

1.  A  brief  quantitative  experimental  investigation  of  the 
effect  of  tilting  the  uniform  and  nonuniform  jets  from  the 
normal  impingement  condition  should  be  made. 

2.  Experimental  measurements  of  the  aerodynamic  forces 

ou  objects  in  grounc’. plane  flows  similar  to  those  observed 
under  the  uniform  and  nonuniform  jets  should  be  made. 

The  resulting  data  should  then  be  combined  with  >,he  com¬ 
pleted  studies  of  the  flow  field  under  impinging  jets  to 
provide  a  means  of  predicting  particle  entrainment  under 
hovering  VTOL  vehicles. 
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TECHNICAL  DISCUSSION 
I  INTRODUCTION 


A  major  attraction  of  VTOL  aircraft  is  their  ability  to  operate  from 
unprepared  sites.  Such  an  ability  does  not  coro'^  without  some  serious 
problems,  one  of  which  is  the  entrainment  of  particles  on  the  ground 
into  the  jet  or  rotor  downwash  during  takeoff,  and  hovering  or  slow- 
speed  flight  near  the  ground.  In  order  to  obtain  an  understanding  of 
the  ground  particle  entrainment  process,  it  is  necessary  to  know  both 
the  nature  of  the  flow  ne.'r  the  ground  and  the  forces  which  act  on 
particles  immersed  in  such  a  flow.  The  research  described  in  this 
report  is  part  of  a  continuing  program  conducted  at  Cornell  Aeronautic 
cal  Laboratory,  designed  to  investigate  particle  entrainment  under 
hovering  VTOL  vehicles. 

To  date,  the  properties  of  the  flow  under  normally  impinging  uniform 
and  nonuniform  jets  have  been  studied.  The  experimeital  investigation 
of  the  flow  in  a  uniform  impinging  jet  and  the  formulation  of  a  method 
of  solution  for  the  equivalent  inviscid  problem  have  been  presented 
in  References  1  and  2.  This  report  presents  the  results  of  an  experi* 
mental  study  of  a  normally  impinging  nonuniform  jet  and  the  completed 
inviscid  analysis  of  the  uniform  impinging  jet.  The  nonuniform  jet  had 
an  approximately  triangular  velocity  profile  similar  to  that  generated 
by  rotors  and  ducted  fans.  In  addition,  an  approximate  method  of 
analysis  to  determine  the  properties  of  the  flow  along  the  ground  after 
the  jet  has  been  turned  is  presented  and  the  results  are  compared  with 
experiment. 
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II  EXPERIMENTAL  PROGRAM 


A.  Experimental  Apparatus 

The  CAL-Air  Force  One-Fbot  High-Speed  Wind  Tunnel  was  used 
during  the  experimental  research.  The  transonic -supersonic  test 
section  was  completely  removed.  A  1 2 -inch-diameter  molded  fiber¬ 
glass  nozzle,  designed  to  produce  a  uniform  parallel  flow  at  the  exit, 
was  mounted  to  an  adapter  plate  on  the  test  section  face  of  the  wind 
tunnel  settling  chamber.  An  axisymmetric  shear  screen  mounted  in  a 
circular  1 -foot-diameter  tube  (Figure  1)  was  fastened  to  the  end  of 
the  uniform  jet  nozzle.  This  shear  screen  was  designed  to  generate 
a  near  triangular  jet  velocity  distribution  across  the  radius  at  the  exit 
of  the  extended  nozzle.  The  ground  board  was  8  by  8  feet  square 
and  could  be  set  at  varying  distances  from  the  nozzle.  A  pattern  of 
static  pressure  taps  was  located  on  the  ground  board.  These  taps 
were  connected  to  an  inclined  manometer  board.  Photographs  of 
the  experimental  apparatus,  excluding  the  shear  screen  assembly,  have 
been  presented  in  References  1  and  2. 

Velocities  and  static  pressures  in  the  jet  were  measured  with  total 
and  static  pressure  probes  mounted  on  a  traverse  mechanism.  In 
addition,  some  velocity  profiles  in  the  jet  near  the  ground  were  mea¬ 
sured  with  a  hot-wire  anemometer.  Ground  flow  measurements  were 
obtained  with  small  static  and  total-head  probes  extending  through  holes 
in  the  ground  board  and  mounted  to  a  second  traversing  mechanism,  so 
made  tha'  the  position  of  the  probe  relati  'o  to  the  ground  could  be 
determined  with  accuracy.  The  boundary  layer  total-head  probe  was 
made  from  0.014-inch  O.  D.  hypodermic  tubing  with  a  tip  tapered  to 
0.  009  inch.  The  static  probe  was  made  from  hypodermic  tuning  with  an 
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outside  diameter  of  0.  031^  inch.  The  boundary  layer  probes  were 
attached  to  a  small  U-tube  manometer  on  the  reverse  side  of  the 
ground  board. 

A  flow  visualization  technique  was  used  to  obtain  pictures  of  the  flow 
in  the  impinging  jets.  A  typical  arrangement  for  these  tests  is  shown 
in  Figure  2.  This  photograph  was  taken  during  the  earlier  un  'orm 
jet  experiments, so  the  shear  screen  is  not  shown.  In  these  tests,  a 
thin  steel  or  aluminum  plate  (splitter  plate)  was  suspended  by  an  angle* 
iron  frame  in  the  horizontal  plane  of  symmetry  of  the  jet.  The  frame 
was  constructed  so  that  it  did  not  interfere  with  the  flow.  The  plate 
was  coated  with  a  mixture  of  powdered  graphite  (lampblack)  and 
kerosene, and  the  flow  was  started.  When  a  well-defined  pattern  had 
developed, the  flow  was  stopped.  When  the  mixture  was  completely  dry. 
the  plate  was  removed  and  photographed.  In  addition,  patterns  were 
obtained  for  the  flow  along  the  ground  under  the  nonuniform  jet.  No 
flow  visualization  pictures  were  recorded  for  the  flow  along  the  ground 
under  the  uniform  jet  because  exploratory  tests  showed  nothing  of 
significance,  only  the  expected  symmetrical  pattern  of  radial  stream¬ 
lines  originating  at  the  jet  centerline. 

B.  Experimental  Results 

1 .  Flow  Visualization  Pictures 

Since  the  lampblack  and  kerosene  flow  pictures  provide  a  visualization 
of  the  turning  process, these  pictures  will  be  presented  first  so  that  the 
quantitative  measurements  can  more  easily  be  interpreted.  Some 
general  comments  on  the  technique  and  its  limitations  are  given  in 
Appendix  1.  A  previously  obtained  flow  picture  for  the  uniform  im¬ 
pinging  jet  is  included  for  purposes  of  comparison. 
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Photographs  of  the  patterns  obtained  are  presented  in  Figures  3  through 
5  for  H/D  =  2.  The  streamline  pattern  on  the  splitter  plate  in  the  uniform 
jet  (Figure  3)  behaves  as  expected.  Near  the  nozzle,  there  is  a  central  area 
of  potential  flow  surrounded  by  a  mixing  region  which  grows  in  thickness 
with  distance  from  the  nozzle.  Turning  of  the  jet  appears  to  occur  first 
at  a  distance  of  approximately  one  nozzle  diameter  above  the  ground.  In 
this  picture,  the  trailing  edge  of  the  plate  was  6  inches  from  the  ground 
level.  This  rather  large  gap  along  with  the  cutout  in  the  rear  center  por¬ 
tion  of  the  splitter  plate  was  required  to  prevent  boundary  layer  separation. 
This  effect  is  discussed  in  Appendix  I. 

With  the  nonuniform  jet  the  flow  in  the  deflecting  portion  of  the  jet  is  altered 
considerably  (Figure  4(a)  ).  The  trailing  edge  of  the  splitter  plate  in  this 
picture  is  only  1  inch  from  the  ground.  There  is  a  dome-shaped  region 
surrounding  the  center  of  the  impingement  area  in  which  the  flow  moves 
radially  inward  along  the  ground  and  then  curves  upward  and  outward 
to  mix  with  the  impinging  flow.  At  the  base  of  this  dome, there  is  a  toroidal- 
or  doughnut  -shaped  vortex  sitting  a  short  distance  above  the  ground.  At 
large  radial  distances,  the  flow  near  the  ground  approached  that  of  the 
classical  turbulent  radial-wall  jet,  just  as  in  the  case  of  the  uniform  jet. 

The  flow  under  the  nonuniform  jet  is  similar  to  that  found  under  annular - 
jet  ground-effect  machines.  Note  that  there  is  a  region  of  upwash  under 
jets  having  a  velocity  which  increases  with  the  distance  from  the  axis. 

The  flow  along  the  ground  which  corresponds  to  Figure  4  (a)  is  shown 
in  Figure  4  (b).  The  inflow  region  is  surrounded  by  a  stagnatio;: 
ring  (  dark  circle  )  where  the  velocity  is  essentially  zero.  Outside  of 

*  The  black  square  in  the  center  of  the  photograph  is  a  piece  of  blotting 
paper  which  was  used  to  soak  up  an  excess  oi  lampblack  and  kerosene 
mixture  which  tended  to  accumulate  in  that  region.  Deviation  of  the 
indicated  streamline  pattern  from  radial  flow  just  above  the  blotting 
paper  is  a  gravity  effect  on  the  kerosene-lampblack.  This  occurred 
because  the  board  was  mounted  vertically  in  the  tests. 
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this  ring,  the  flow  is  directed  rauially  outward.  In  all  cases  for  which 
flow  pictures  were  obtained  (H/O  -  4,  2,  1,  the  stagnation  ring 
coincided  with  the  position  of  maxinnum  static  pressure  along  the  ground, 
and  its  radius  was  larger  than  the  radius  of  the  toroidal  vortex  observed 
on  the  corresponding  splitter -plate  pictures. 

The  general  features  of  the  pictures  obtained  at  the  other  values  of 
H/D  were  similar  to  those  shown  here,  the  only  difference  being  that  the 
size  of  the  dome- shaped  region,  the  toroidal  vortex  and  the  stagnation  ring 

all  became  smaller  as  H/D  was  either  increased  or  decreased  from  the 
value  of  2.  A  probable  explanation  for  this  effect  is  chat  on  one  hand, 
as  H/D  is  decreased,  the  jet  is  constrained  by  the  nozzle  walls  and 
so  the  region  in  which  turning  of  the  jet  occurs  must  be  reduced.  On 
the  other  hand,  as  H/D  is  increased  beyond  2,  mixing  of  the  jet  with 
still  air  prior  to  turning  tends  to  reduce  the  effective  shear  in  the 
jet  so  that  any  effect  associated  with  shear  in  the  jet  before  impinge¬ 
ment  will  also  tend  to  be  reduced.  The  type  of  shear  being  considered 
here  is  that  which  will  provide  a  velocity  minimum  on  the  jet  axis. 
Parabolic  profiles  such  as  are  found  iii  fully-developed  pipe  flow  have 
a  shear  of  the  opposite  sense.  This  latter  type  of  velocity  profile  will 
not  form  a  toroidal  vortex  with  its  associated  central  region  of  reversed 
flow. 

Figure  b  illustrates  the  effect  of  small  deviations  from  the  normal 
impingement  condition  for  the  nonuniform  jet.  These  photographs  were 
obtained  with  a  tilt  angle  of  5  degrees  between  the  nozzle  axis  and  normal 
to  the  ground  board.  The  distance  between  the  nozzle  and  ground  in  this 
case  was  measured  along  the  extended  axis  of  the  nozzle.  The  distance 
between  the  trailing  edge  of  the  splitter  plate  and  the  ground  was  1  inch. 
The  white  spot  in  Figure  5(b)  resulted  from  wiping  off  an  excess  of 
lampblack  and  kerosene  which  accumulated  at  that  point.  Kven  for  this 
relatively  small  angle  of  tilt,  the  flow  pattern  has  changed  considerably 
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from  that  obtained  with  the  nonuniform  jet  impinging  normally.  The 
flow  near  the  ground  is  appreciably  thicker  on  the  side  furthest  removed 
from  the  nozzle  (right  side  cf  Figure  5(a)  ),  and  the  stagnation  ring  has 
been  broken  on  this  same  side  (left  side  of  Figure  5(b)  ).  Although 
no  quantitative  measurements  were  made  with  the  tilted  jet,  it  is 
apparent  from  the  pictures  alone  that  the  nonuniform  jet  is  extremely 
sensitive  to  tilt  angle.  No  check  of  the  effect  of  tilt  angle  was  made 
during  the  earlier  uniform  jet  experiments,  but  it  is  believed,  on  the 
basis  of  the  general  experience  gained  in  setting  up  conditions  for 
normal  flow  impingement  for  both  types  of  jets,  that  the  uniform  jet 
is  not  as  sensitive  as  the  nonuniform  jet  to  small  angles  of  tilt. 

2.  Quantitative  Measurements 

Tests  were  conducted  on  the  nonuniform  jet  at  ground-board  locations 
corresponding  to  nozzle-height-to-nozzleodiameter  ratios,  H/D,  of 
4.  07,  1. 99  and  0.  49.  Most  of  the  data  were  obtained  with  the  settling 

chamber  pressure  held  constant  at  16  inches  of  water,  which  provided  a 
maximum  velocity  near  ♦he  outer  edge  of  the  jet  of  approximately 
188  fps  at  the  nozzle  exit.  Sample  measurements  were  also  made 
with  a  maximum  jet  velocity  of  119  fps  at  the  nozzle  exit.  At  each 
H/D,  velocity  surveys  were  made  at  the  nozzle  exit  and  in  the  flow 
along  the  ground  at  various  radial  distances  from  the  stagnation  point 
up  to  a  maximum  of  r  =  42  inches.  Static  pressure  distributions  were 
measured  along  the  jet  centerline  between  the  ground  board  and  the  jet 
nozzle  and  also  along  the  ground  board.  The  notation  used  in  the  pre¬ 
sentation  of  the  data  is  shown  in  Figure  6. 

It  was  found  that  the  flow  at  the  nozzle-exit  plane  was  somewhat 
irregular.  A  typical  set  of  velocity  data  is  shown  in  Figure  7.  This 
figure  shows  the  velocity  measured  along  four  radial  lines  at  90-degree 


intervals  in  the  exit  plane  of  the  nozzle  for  the  jet  exhausting  into 
essentially  free  air.  The  data  have  been  made  nondimensional  by  a 
reference  velocity  44f»which  was  chosen  for  experimental  convenience 
and  which  is  a  reflection  of  the  pressure  in  the  settling  chamber  of 
the  wind  tunnel.  The  dashed  line  is  an  estimated  mean  curve  drawn 
through  the  data.  This  curve  was  used  to  evaluate  the  momentum  and 
mass  flow  of  the  jet.  The  accuracy  of  such  a  procedure  can  be 
judged  by  comparison  of  the  jet  momentum  calculated  from  mean 
curves  for  different  distances  from  the  nozzle.  The  nondimensional 
jet  momentum.  •  calculated  from  the  mean  profile  at  the 

nozzle  exit,  was  0.  302,  as  compared  to  0.  308  for  =  2.  02.  This 
is  a  difference  of  2  percent,  which  is  quite  satisfactory. 

It  can  be  seen  that  the  triangular  velocity  piofile  for  which  the 
shear  screen  was  designed  was  only  partially  attained.  In  particular, 
the  velocity  at  the  outer  edge  of  the  jet  deviates  considerably  from 
the  desired  distribution.  The  irregularities  in  this  region  arc  a 
consequence  of  insufficient  nozzle  length.  A  greater  length  of 
nozzle  downstream  of  the  shear  screen  would  have  allowed  more 
mixing  of  the  flow  and,  hence,  smoothed  out  the  velocity  profile. 

The  greater  nozzle  length  was  prohibited,  however,  by  the  space 
available  for  the  tests. 

The  mean-velocity  profiles  in  tl.e  free  nonuniform  jet  measured  at 
various  distances  from  the  nozzle  arc  shown  in  Figure  8.  The 
profile  of  a  uniform  jet  with  equal  thrust  is  shown  for  comparison. 

The  irregularities  in  the  profiles  disappeared  rapidly  with  distance 
from  the  nozzle.  At  the  same  time  the  peak  velocity  in  the  jet 
moved  toward  the  axis  and  decayed  almost  linearly  with  distance 
from  the  nozzle,  while  the  velocity  along  the  jet  centerline  increased. 
The  behavior  of  the  maximum  velocity  in  the  free  jet  and  the  velocity 
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on  the  jet  centerline  are  shown  in  Figure  9.  The  decay  of  the  maxi> 
mum  velocity  for  this  particular  nonuniform  jet  is  represented  quite 
well  by  the  expression  =  0.  695  -  0.0225  •*  . 


A  comparison  of  the  entrainment  in  the  uniform  and  nonuniform  jets 
is  shown  in  Figure  10.  In  this  figure  the  mass  flow  in  each  jet  is 
plotted  as  a  function  of  distance  from  the  nozzle.  The  maos  flow  has 
been  normalised  by  the  value  measured  at  the  nozzle  exit.  When 
reduced  in  this  way,  the  entrainment  properties  of  the  two  jets  appear 
to  be  identical.  The  dashed  line  is  an  empirical  straight-line 
approximation  to  the  data. 


Figure  11  shows  ti  e  velocity  distributions  at  the  nozzle  for  H/D  =  4, 

2,  1  and  0.  5;  the  velocities  are  again  normalized  by  the  reference 
velocity  .  The  velocity  distribution  for  H/D  =  4  was  essentially 
identical  to  that  for  H/0  =  «•  .  The  effect  of  the  proximity  of  the 
ground  is  most  evident  near  the  jet  axis.  Perhaps  a  better  indication 
of  the  ground  effect  is  given  by  Figure  12,  where  the  static -pressure 
profiles  at  the  nozzle  are  plotted.  The  pressures  have  been  nondimen- 
sionalized  by  the  settling  chamber  pressure,  .  The  static  pressure 
difference  was  zero  for  H/D  =  4  and  ••  .  At  H/D  =  0.  49,  the  non- 
dimensional  static  pressure  on  the  jet  axis  reached  a  value  of  0.25. 
This  is  approximately  50  percent  of  the  maximum  total  pressure 
measured  at  the  outer  edge  of  the  jet.  There  was  some  ground  effect 
even  at  H/D  =  1. 99.  The  corresponding  tests  c(  the  uniform  jet 
showed  no  back  pressure  at  H/  D  =  2.  The  first  evidence  of  ground 
effect  in  that  case  appeared  at  H/D  =  1. 

The  static  pressure  distributions  on  the  ground  are  presented  in 
Figures  13(a)  through  13(d).  Data  are  shown  for  four  radial  lines 
separated  by  90-degree  intervals.  The  line  &  -  270  degrees^  i*  the 

See  List  of  Symbols 
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one  on  which  boundary -layer  meaturementa  were  made.  It  is 
evident  from  the  data  that  perfect  axial  synunetry  was  not  attained 
in  these  experiments.  There  is  a  general  tendency  for  the  peak 
static  pressure  to  be  higher  in  the  horisontal  plane  (  ^  =  90  and 
270  degrees)  than  in  the  vertical  plane  (tf  -  0  and  180  degress). 

The  position  of  the  maximum  static  pressure  coincided  with  the 
position  of  the  stagnation  ring  in  the  flow  visualisation  pictures  as 
one  might  expect.  The  static  pressure  on  the  ground  under  the 
tilted  nonuniform  jet  is  shown  in  Figure  14.  The  plane  correspond¬ 
ing  to  ^  -  90  and  270  degrees  is  the  one  in  which  a  5 -degree 
angle  of  tilt  was  imposed,  with  the  line  B  =  270  degrees  being 
farthest  from  the  nozsle.  The  break  in  the  stagnation  ring  illustrated 
in  Figure  5(b)  is  reflected  in  Figure  14  as  a  suppression  of  the 
peak  in  static  pressure  along  the  line  B  •  270  degrees. 

Figure  15  showi  the  static  pressure  distribution  along  the  Jet 
centerline  as  a  function  of  nondimens ional  distance  from  the 
ground,  */D.  Inspection  of  Figures  13  and  15  shows  that  ihe  static 
pressure  is  approximately  constant  in  the  region  surrounding  the 
center  of  impingement,  and,  in  addition,  this  constant  pressure 
region  is  largest  in  extent  for  H/0  «  2.  These  observations  are 
in  accord  with  the  behavior  of  the  jet  as  determined  from  the  flow 
visualization  pictures.  The  increase  in  static  pressure  along  the 
jet  axis  near  the  ground  at  H/D  =  4  over  that  at  H/O  s  2  is  a  reflec¬ 
tion  of  the  mixing  process  before  impingement.  The  dynamic 
pressure  on  the  centerline  of  the  free  jet  increases  with  distance 
from  the  nozzle  because  of  viscous  mixing  (Figure  8),  and  this 
increase  causes  a  rise  in  static  pressure  recovery  near  the  center 
of  impingement. 


A  brief  hot-wire  survey  of  the  impingement  region  was  performed 
for  the  normally  impinging  nonuniform  jet  with  H/D  =  1.99.  The 
result  of  this  survey  is  shown  in  Figure  16.  The  velocities  shown 
are  absolute  velocities;  the  direction  of  the  flow  is  not  indicated  by 
a  single  hot  wire.  This  should  be  kept  in  mind  when  studying 
Figure  16.  In  particular,  the  flow  visualization  pictures  indicated 
that  the  direction  of  the  velocity  along  the  jet  axis  near  the  board 
was  perpendicular  to  and  away  from  the  board.  The  main  point 
which  can  be  determined  from  Figure  16  is  that  the  velocity  in  the 
dome-shaped  region  of  Figure  5(a)  is  only  a  small  fraction  of  that 
found  outside  the  dome  and,  hence,  should  not  play  an  important 
role  in  the  entrainment  problem.  There  is,  however,  the  possibil¬ 
ity  that  particles  entrained  in  the  flow'  in  higher  velocity  regions 
may  be  carried  to  the  center  of  the  impingement  region  by  this 
secondary  flow  and  deposited  as  a  mound  near  the  axis  of  symmetry 
of  the  flow. 

Figure  17  presents  velocity  profiles  in  the  flow  along  the  ground 
at  various  radial  stations  in  the  turning  region  for  H/D  -  0.49, 

1. 99,  and  4.07.  A  velocity  profile  measured  inside  the  stagnation 
ring  with  H/D  =  1. 99  is  included  in  Figure  17(b).  Note  that  for  radii 
near  the  stagnation  ring,  the  maximum  velocity  is  reached  at  a 
considerable  height  above  the  ground.  Below  this  maximum  velocity 
point,  the  velocity  decreases  gradually  as  the  wall  is  approached 
until  at  a  very  small  distance  from  the  wall,  the  beginning  of  a 
true  boundary  layer  is  discernable.  The  region  of  increasing 
velocity  between  the  very  thin  boundary  layer  and  the  maximum 
velocity  is  due  to  the  static  pressure  difference  between  the  wall 
and  the  outer  edge  of  the  flow.  The  decay  in  veloc  ity  at  larger 
distances  from  the  wall  is  due  to  viscous  mixing.  As  the  radius  is 
increased  to  regions  w'here  the  static  pressure  gradient  at  the  wall 
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approaches  zero  (see  Figure  13),  the  velocity  profiles  rapidly  approach 
the  fo.  m  associated  with  the  classical  turbulent  radial  wall  jet. 

The  experimental  data  tor  H/D  =  1.99  and  r/D  >  1  are  plotted  as  u/U^^ 
vs  y/y^  in  Figure  18;  is  the  distance  to  the  point  in  the  outer 
mixing  region  at  which  u,  =  Figure  6).  In  general,  the  data 

fall  on  a  single  curve  when  plotted  in  this  way  so  long  as  the  radius  is 
sufficiently  large  that  the  pressure  gradient  at  the  ground  is  approxi¬ 
mately  zero.  Identical  results  were  obtained  for  H/D  -  0.49  and  4.  07 
with  the  exception  that  the  minimum  radius  at  which  the  data  collapsed 
to  a  wall-jet  type  of  velocity  profile  changed  with  H/D.  For  each  H/D, 
the  largest  value  of  r/D  indicated  in  Figure  17  is  the  smallest  r/D  at 
which  the  wall -jet  profile  was  obtained.  The  next  smaller  value  of  r/D 
in  each  case  followed  the  wall-jet  profile  in  the  outer  portion  of  the  flow 
except  for  H/D  ~  0.49  ( see, e.  g. ,  Figure  18,  r/D  s  1.  33).  Sample  data 
at  a  second,  lower,  value  of  mass  flow  in  each  case  gave  essentially 
the  same  results  when  nondimensionalized.  This  is  illustrated  in  Figure 
18  for  r/D  =  2.10. 

«* 

In  Figure  19  the  boundary  layer  region  for  the  data  of  Figure  18  is 
plotted  on  a  greatly  expanded  scale.  The  cc  mments  made  above  apply 
to  Figure  19  as  well. 
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m  THEORETICAL  ANALYSES 


A.  Uniform  Inviscid  Impinging  Jet  Theory 

A  theory  for  the  flow  along  the  ground  plane  of  an  impinging  jet,  particu¬ 
larly  for  the  boundary  layer  flow,  requires,  first,  that  the  ground-plane 
pressure  distribution  be  predicted.  It  was  recognized  that  for  free  jets 
in  air,  viscous  effects  are  important.  However,  the  work  reported  in 
Reference  1  indicated  that,  for  H/D  less  than  approximately  2,  the 
ground-plane  pressure  distribution  under  the  uniform  impinging  jet  was 
essentially  that  of  the  inviscid  jet.  Furthermore,  the  available  theoreti¬ 
cal  results  for  inviscid  impinging  jets  were  not  applicable  for  H/D  less 
than  2.  For  H/D  less  than  2,  as  the  jet  exit  approaches  the  ground  plane, 
there  is  a  back-pressure  which  is  reflected  at  the  exit  in  terms  of  re¬ 
duced  dynamic  pressure  in  the  center  of  the  jet.  Hence,  the  jet  exit 
velocity  profile  is  no  longer  uniform.  As  H/D  is  reduced  below  2,  there 
are,  in  addition,  substantial  changes  in  t..  e  ground-plane  pressure  distri¬ 
bution. 

analysis  of  the  uniform  inviscid  impinging  jet  was  a  continuation 
of  me  study  partially  reported  in  Reference  1.  Work  was  brought  to  a 
conclusion  during  the  present  research.  The  analysis  is  based  on  a 
vortex  sheet  representation  for  the  jet  boundaries.  The  resultant 
mathematical  model  was  programmed  for  an  IBM  704  digital  computer. 

An  iterative  technique  was  developed  in  which  successive  positions  of  the 
jet  boundary  were  computed  until  the  appropriate  boundary  conditions 
for  the  flow  were  satisfied.  Details  of  the  mathematical  model  and  the 
implementation  of  this  model  on  the  computer  are  presented  in  Appendix  II. 
For  the  sake  of  completeness,  parts  of  the  discussion  of  Reference  1  are 
repeated. 
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During  the  formulation  stage  of  the  analytical  program,  it  was 
r;.alized  that  problems  could  arise  concerning  convergence  of  the 
iterative  technique  which  might  be  amplified  by  the  approximations 
inherent  in  any  digital  computer  program  for  continuous  systems. 
Hence,  the  convergence  properties  were  investigated  from  an  experi¬ 
mental  point  of  view  by  a  trial  and  error  process.  Fortunately,  con¬ 
vergence  was  ultimately  attained. 

Boundaries  and  flow  properties  were  computed  for  H/D  =  1  and 
H/D  =  0.  25.  The  computed  boundaries  are  shown  in  Figure  20 . 
where  the  results  of  Leclerc,  Reference  3  ,  for  H/0  =  2  (approxi* 
mately  equivalent  to  H/D  s  ss,are  shown  for  comparison.  Computed 
static  pressure  distributions  along  the  jet  centerline  are  shown  in 
Figure  21  ,  and  along  the  ground  plane  in  Figure  22  .  Computed 
velocity  profiles  and  flow  angularity  with  reference  to  the  jet  axis 
at  the  nozzle  are  presented  in  Figure  23. 

B.  Boundary  Layer  Analysis 


An  approximate  method  of  analysis  was  developed  in  Reference  1 
which  appeared  to  he  adequate  for  predicting  the  flow  along  the 
ground  under  an  impinging  uniform  jet.  The  method  was  semi- 
cmpirical  in  that  it  was  based  on  the  experimentally  determined  fact 
that  the  mass  flow  along  the  ground  could  be  made  independent  of 
H/D  (  at  least  in  the  range  H/D 2 4)  by  a  proper  choice  of  non- 
dimensionalizing  parameters.  In  addition,  the  analysis  required 
knowledge  of  the  pressure  distribution  along  the  ground.  This  latter 
requirement  was  to  be  met  by  the  inviscid  solution  to  the  problem  of 
the  impinging  uniform  jet.  Such  a  solution  is  reported  in  Section  Ill  (A). 
Since  this  solution  was  not  available  at  the  time  the  boundary  layer 
analysis  was  performed,  the  experimental  pressure  distributions  were 
used  in  anticipation  of  the  invisiid  solution. 
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1  •  ,  .  f  I  t  '  t  ^  ‘  *  ’  .  ’  :  '  , j  ;  ‘  * 

The  problem  became  more  difficult  when  an  initially  nonupiform,  jet 

j  ,  ;  (  '  ‘  .  I  .  ,  1  .  I  '  ’  '  j  ^  ‘  ‘  '  '  * 

was  considered.  Since  a  theoretical  solution  for  the  pre*  sure  distri* 

>  :■  ,  ‘if,..:-  M-  j  .  ■  .  '  • 

bution  on  the  ground  under  such  a  jet  was  not  available,  it  was 

I . ' !  ■  r . ) ,  1 *1  ,  1  '  '  I  j  ■  j .  •  '  ■  ■  ■  '  • 

necessary  to  restrict  the  analysis  to  regions  of  the  flow  where  the 

.•  I .  'T  ’  V  '  H  f  "  >  C|  .i  ^  ’  ;  ■ '  • ;  J  ■  1  I 

pressure  gradient  is  approximately  zero.  The  logical  proceuure 

:  . .  !  •  II  :  >  !  . .  *  ■  I ’  V.'  f  ^  :  1  i  •  .  1 1  ;  I  •  . . 

appeared  to  be  to  search  for  a  correlation  between  the  mass  flow  param* 
cters  of  the  uniform  and  nonuniform  jets  and  to  use  such  a  correlation, 
if  found,  in  conjunction  with  a  simple  momentum  balance  to  describe 
the  flow  at  radii  sufficiently  large  for  the  pressure  to  have  reachet^. 

I  r  *  ^  ‘ 

nearly  ambient  conditions^  A  metl^od  of  analysis  based  on  this 
correlation  would  then  be  applicable  to  both  the  uniform  and  the  pon- 
uniform  jets. 


.  j  '  ' 

A  factor  in  favor  of  this  procedure  is  that  the  velocity  profiles  near  the 
ground  in  the  nearly  ambient  pressure  region  can  be  found  from  Glauert's 
radial  wall  jet  analysis  (Reference  4  i.  The  turbulent  flov'  solution  is 
the  one  of  interest  for  this  investigation.  The  shape  of  the  theoretical 

.  I  ; 

velocity  profile  was  prescribed  by  a  parameters,  which  in  turn 

ti  s 

depended  weakly  on  the  Reynold's  number  where  is  the  Ipcc'  1 

peak  velocity  in  the  radial  flow  and  6^  is  the  distance  between  the 
point  at  which  u  s  and  the  point  at  which  as/t^lFigure  6). 


It  is  desirable  to  obtain  an  estimate  of  Ct  from  the  properties  of  the 
jet  before  impingement.  Because  of  the  weak  dependence  of  on 

,  this  can  be  accomplished  with  considerable  accuracy  ever,  if 

y 

some  relatively  broad  assumptions  are  made,  it  was  assumed  that 
Glauert's  theory  applies  for  r/R  -  3,  that  the  loss  in  radial  momentum 
flux  due  to  skin  friction  is  negligible  at  this  radius,  and  that  equals 
the  value  of  in  the  free  jet  at  a  distance  t/R*RIR  f )  where  •^•3 

for  this  calculation.  A  discussion  of  this  method  of  estimating  cc  has 

1 

been  presented  in  Appendix  II  of  Reference  I  and  will  not  be  repeated 
here.  The  value  dS«l.  16  computed  in  this  manner  compares  very  tlosely 
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with  those  derived  from  the  experimental  measurements  of 


Once  the  shape  of  the  velocity  profile  has  been  determined  by  Ob  ,  the 
functional  dependence  of  and  on  r  is  given  by 


^m/^m  -  A  y,/s  iriR) 


f  J  *  l\. 


(i) 


^  where  the  exponents  "a"  and  "b"  depend  on  the  value  of  <t ,  the  constants 
and  may  be  functions,  of  H/D  and  the  type  of  jet  considered, 
and  CId  is  a  reference  velocity  which  will  be  defined  shortly.  The  approxi 
mate  detern\ination  of  an<t  ^  forms  the  remaining  portion  of  the 
analysis. 


Consideration  of  the  momentum  in  the  jet  provides  one  relation  for 
finding  the  unknown  constants  and  (y .  Consider  an  impinging  jet  in 
which  both  the  total  and  the  static  pressure  distributions  at  the  nossle 
exit  are  functions  of  radius  (Figure  24).  lf,for  the  moment,  one  arti¬ 
ficially  allows  no  mixing  to  occur,  then  one  may  Say  that  the  total 
pressure  along  any  streamline,  will  be  a  constant.'  Along  any  strei  m 
tube  annulus  (dotted  lint;  in  Figure  24)  in  the  impinging  >et. 


Zrryo  *  Zn/o  r  u(^)dff , 

Pr  “  ~  constant  in  the  streamtube  and 


I  V  ‘  '  1 


(2) 


(3) 


where  the  subscript  .>V  denotes. that  the  quantities  are  evaluated  at  Jtbe 
nozzle  exit  plane.  If  one  considers  the  flow  along  the  ground  where 
Pm  -P*0,  then 


a(^)  .  li 


/> 


IZ<Pt  niM 

i-y—' 


(4) 


T^herefOre  from  (2),  (3  ),  and  (4) ,  * 

-771755 - T  y  -(Pt-PaU - 
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or 


an^ 


u 


¥ 


.m 


(5) 


Note  that  the  thickness  of  the  flow,  y'  (r)  $  is  obtained  simply  by  replac- 

s 


P^H  .  t* 
mg  /  with 


in  equation  (5).  The  radial  momentum  flux  at  r  is 


given  by 


*  ZirprJ  u‘dy. 


(6) 


Substitution  of  equation  (5)  into  (6)  gives 


(7) 


which  is  the  radial  momentum  flux  in  terms  of  conditions  at  the  nozzle 
exit  plane  provided  that  the  radial  station  considered  is  sufficiently  far 
from  the  jet  axis  so  that  the  pressure  at  the  wall  has  become  approxi¬ 
mately  atmospheric. 


Since  Equation  (7)  deals  with  momentum  flux,  it  is  valid  for  the  flow 
under  a  real  impinging  jet  even  though  entrainment  was  neglected  in 
its  derivation,  provided  that  the  momentum  flux  loss  due  to  wall  skin 
friction  can  be  neglected.  The  momentum  flux  loss  at  the  start  of  the 
wall-jet  flow  has  been  found  experimentally  to  be  small  for  all  cases 
tested  to  date,  so  Equation  (7)  will  be  used  in  the  following  analysis. 


It  is  convenient  at  this  point  to  define  two  characteristic  mean  velocities 
based  on  conditions  at  the  nozzle  exit.  These  are  a  mean  mass  flow 
velocity 


nr 

and  a  mean  momentum  velocity  ^  | 

With  this  definition  of  IJ^  .  Equation  (7) 


where  i/; 


t  •  m  m  m  m mm 

.  -  JIK 

"  I 
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simplifies  to 


_ ^ 

.e  ■>  ~  f 


(skin  friction  neglected) 


(7a) 


An  additional  relationship  is  required  to  obtain  the  two  unknowns, 

^  and  £2:)uation  (1).  In  the  tests  of  the  uniform  jet  (Reference 

1  ).  this  relation  was  obtained  empirically.  It  was  found  that  the 

radial  mass  flow,  •  in  the  wall- iet  region  was  independent  of 

Or  /U 

H/D  when  plotted  as  jrppiujCWt'^R  where  is  the  velocity  at 

the  edge  of  the  uniform  jet.  On  reviewing  the  calculations  for  the 
uniform  jet,  however,  a  small  error  in  the  calculation  of  u/u^  from 
the  experimental  data  was  found  for  H/D  =  0.  25.  The  revised 
estimate  of  became  0.  56  as  compared  to  the  original  estimate 

of  0.  60.  With  this  more  accurate  value  for  jf-  ,  the  nondimensional 

Qf.  /c/ 

mass  flow  parameter  became  \~^l'  revised  parameter 

will  be  used  in  the  following  analysis.  Application  of  the  definition 
of  to  the  uniform  jet  (  note  *(4  the  uniform  jet)  shows  that 
that  the  equivalent  mass  flov;  parameter  in  terms  of 

Qft  /t7 

Uta  is  correlation  of  the  experimental  data 

for  both  the  uniform  and  the  nonuniform  jets  is  illustrated  in  Figure  25. 


The  data  for  each  type  of  jet  correlate  when  plotted  in  this  way  but  there 
is  15  percent  difference  between  the  two  sets  of  data.  It  is  possible 
that  some  of  the  difference  may  be  attributable  to  small  deviations 
from  axial  symmetry  in  the  te<ts  of  the  nonuniform  jet.  Such  devi¬ 
ations  are  noted  in  Sr^ction  11.  A  further  indication  of  asymmetry  in 
the  flow  under  the  nonuniform  jet  is  that  the  experimental  radial 
momentum  flux  with  H/D  =  4  reached  values  up  to  10  peicent  higher  than 
is  theoretically  possible  for  an  axisymmetric  jet  with  no  losses. 

Since  the*  measurements  in  the  flow  along  the  ground  were  made  only 
along  one  radial  line,  the  results  are  subject  to  error  if  the  flow  is 
not  accurately  axisymmetric.  In  view  of  this  uncertainty  in  the 
experimental  results,  it  was  decided  that  a  weighted  average  of  the 
two  sets  of  data  would  be  used  to  complement  Equation  (7a)  in  the 
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L' a  h  alation  of  C  ,  and  6 


u 


The  variation  of  0,.  with  r  is  prescribed  by  Glauert's  form  param¬ 
eter  cx  .  For  a  =  1.  16,  the  theoretical  values  of  a,  and  b  in 
Equation  (1)  are  -1.075  and  1.008  respectively.  However,  the  ex-  ‘ 

perimental  values  of  a  and  b  as  found  in  the  uniform  jet  tests  proved  • 

to  be  better  represented  by  a  =  -1.  143  and  6=  1. 028.  Since,  in  » 

practice,  the  value  of  a  is  unlikely  to  differ  very  much  from  that 
found  in  the  present  work  because  of  an  insensitivity  of  cc  to  nozzle 
Reynold's  number,  and  the  values  of  a  and  b  in  turn  are  only  weakly  • 

dependent  onczin  this  range,  it  was  felt  justified  to  assume  that  the 
experimental  values  of  ol  and  b  were  representative  of  all  practical 
cases.  Hence,  ~  The  relations 

for  the  uniform  jet 

=  O ^3  for  the  nonuniform  jet 

appear  to  give  the  best  tit  to  the  experimental  data.  A  weighted 

/U  /r\^*** 

average  of  these  two  curves,^— i^J  used  in  the  subse¬ 

quent  ana'ysis.  In  arriving  at  this  value,  most  of  the  emphasis  was 
placed  on  the  uniform  jet  data  because  it  was  believed  that  these  data 
were  the  most  acc  urate. 


The  following  relations  determine  the  flow  near  the  ground: 


- 7=-r  / 

TrpR^U^ 


Or  /U^\Of^  ^  J''\ 

mf)  -‘’Mr) 


r  \0  »»r 


and 


r  \  * 


yt/t 

ir 


f'  \i.oie 


(t) 


-^i<  rt 

R  '  R 

-  ^  Cl 

R  ~  R 


-C  ->  Ll 

R  ~  R 


(7a) 

(8) 

(1) 


where  is  the  radius  where  the  wall-jet  flow  begins.  Further,  it 

)  the 


will  be  assumed  that  at  the  start  of  the  wall  jet,  is  equal  to  the 

t  H 

maximum  velocity  in  an  equivalent  free  jet  at  a  distance  —  s 


\ 
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t 


from  the  nozzle  exit.  For  the  two  jets  tested,  the  above  assumption 
gives 

uniform  jet:  •  f  aM  ^  provided  HfD 

(Jut  “ 

nonuniform  jet:  -  /n  tt.  \\  IT) 

(see  figure  9)  ^  if 

where  is  a  reference  velocity  used  for  experimental  convenience  and 
is  representative  of  the  wind  tunnel  settling  chamber  pressure. 


The  method  of  solution  for  and 
(7(a))  one  obtains 


is  as  follows.  From  Equation 


for  JL  ,  A 
R  R 


and,  since  has  been  assumed  in  this  region,  then 

R  15l\u^)  \RI 


and  ♦ 


or 


Z.22  ^ 


•f 


rr/yR^O^ 


for 

for 


R 

r 


A 

R 

A 

R 


ot^ 


for  _ 


TT/OR' 


r 

R 


The  combination  of  Equations  (8)  and  (11)  gives 

^  /r.  \0-99S 

f. 

I 

which  may  be  solved  for  with  the  aid  of  Equation  (9),  or  its 
equivalent  if  other  jet  profiles  are  of  interest. 


iiai^ud viuiis  \o;  diiti 


(10) 


(11) 


(12) 


Once  q//?  has  been  found,  then  and  are  given  by 


(13) 


The  numerical  constants  in  these  equation  were  obtained  by 
integration  of  the  theoretical  veioctiy  profiles,  That  is, 

0^  =  5:  2.22  rrpru^i^^,  for  1.t<a<  f.2 

and 

»  27rpru^!4y,J^  ^  i  . 
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Equations  (1),  (9),  (10),  (12)  and  (13)  along  with  Glancrt's  velocity 
profile  for  a.  =  1.  16  constitute  a  complete  description  of  the  flow  along 
the  ground  in  the  region  where  the  wall  static  pressure  gradient 
approaches  zero.  The  solution  may  be  applicable  for  radii  slightly  less 
than  .  The  flow  near  the  ground  at  will  have  a  minimum  in 
thickness  and  an  approximate  maximum  in  velocity, and  hence  this  is 
probably  the  radius  where  the  entrainment  problem  will  be  most 
severe. 

The  results  of  the  calculations  are  presented  in  Figures  26  and  27 
for  the  uniform  and  nonuniform  jets  respectively. 
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IV  COMPARISON  OF  THEORY  AND  EXPERIMENT. 
AND  DISCUSSION  OF  RESULTS  ' 


A.  Invitcid  Impinging  Uniform  Jet  Analyis 


Computed  jet  boundaries  are  presented  in  Figure  20  for  H/D  :  1.  0 
and  H/D  =  0.  25.  Also  shown  is  the  boundary  derived  by  Leclerc  from 
an  electric  analog  tank  experiment  for  H/D  =  2.  Comparison  of  these 
computed  boundaries  with  experimental  data  is  not  possible,  because 
the  experimental  jet  boundaries  are  obscured  by  viscous  mixing.  It 
is  apparent  from  Figure  20  that  as  H/D  decreases,  the  theoretical 
thickness  of  the  jet  above  th  2  ground  plane  decreases.  This  was  inferred 
in  Reference  1  from  mass  flow  considerations.  As  H/D  decreases, 
the  equivalent  uniform  velocity,  U  ,  in  the  tube  (or  nozzle)  also  decreases. 
Conservation  of  mass  requires  that  at  some  distance  along  the  ground 
plane  from  the  stagnation  point,  where  the  pressure  has  fallen  to  near 
ambient  pressure. 


v;(r) 


(14) 


which  is  Equation  (1)  of  Reference  1  ;  Y^ir)  is  the  height  of  the  jet 
boundary  above  the  ground  plane.  Equation  (14)  is  consistent  with 
Equation  (5),  Section  111  B,  which  applies  generally  for  impinging  jets 
with  no  mixing. 


Computed  static  pressure  distributions  along  the  jet  centerline  and  along 
the  ground  plane  are  compared  with  corresponding  experimental  results 
from  Reference  1  in  Figures  21and  22  respectively.  It  is  apparent  that 
excellent  agreement  was  obtained,  although  computed  ground  plane 
pressures  are  slightly  low  for  0.  7  <  <  1.  4  at  H/D  =1.0.  The  small 

discrepancies  probably  reflect  the  effects  of  viscosity  (mixing),  in  view 
of  the  excellent  agreement  obtained  for  the  ground  plane  pressure  dis¬ 
tribution  at  H/D  =  0.25,  and  for  the  jet  centerline  pressure  distribuuon 
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for  H/D  =1,0.  Comparison  of  computed  and  experimental  velocity 
profiles  at  the  nozzle  are  shown  in  Figure  23  .  Again,  agreement  is 
good,  particularly  for  H/D  =  0.  25.  Also  shown  in  Figure  23  is  the  com¬ 
puted  variation  of  flow  angularity  at  the  nozzle  with  reference  to  the  jet  axis 
for  H/  D  =  0.25  .  Of  interest  is  the  fact  that  the  computed  angularity  is 
as  large  as  18.  5*in  this  particular  instance. 

Of  those  values  of  H/D  for  which  experimental  data  are  available  from 
Reference  1  ,  the  H/D  =  0.25  case  is  the  one  most  critically  affected 
by  the  ground.  The  excellent  agreement  with  experimental  results  for 
H/D  =  0.2  5  clearly  demonstrates  the  accuracy  of  the  theory. 

Convergence  of  the  iteration  process  was  relatively  slow.  The 
H/D  =  1.0  computation  required  approximately  70  iterations  and  about 
4  hours  of  computing  time  on  the  CAL  IBM  704  computer.  The  iteration 
technique  involves  so-called  "gain"  constants  (See  Equation  (II-5)  of 
Appendix  II).  The  appropriate  gain  constants  for  convergence  must  be 
determined  by  trial  and  error.  If  they  are  too  large,  successive  curves 
computed  as  the  iteration  program  proceeds  tend  to  oscillate  through 
larger  and  larger  amplitudes  wnich  ultimately  diverge  violently.  For 
somewhat  smaller  values,  the  oscillations  become  damped,  and  the 
program  ultimately  converges.  If  the  "gain"  constants  are  reduced 
still  more, the  successively  computed  curves  tend  to  approach  the  solution 
asymptotically.  Apparently,  there  is  an  optimum  "gain"  constant  which 
is  analogous  to  critical  damping  and  for  which  convergence  will  be  obtained 
most  rapidly.  Unfortunately,  a  detailed  investigation  of  this  point  was 
not  believed  to  be  practicable  during  the  present  program.  Once  "gain" 
constants  were  determined  lor  which  the  successive  iterations  showed  a 
converging  solution,  these  constants  were  used  thereafter,  and 

:;c 

the  program  was  permitted  to  run  to  convergence.  It  is  believed  that  the 

An  improved  iteration  technique  which  should  provide'  much  more 

rapid  convergence  is  outlined  in  Appendix  II. 
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convergence  criterion  used  in  these  con^putations  insures  an  accuracy 
of  between  three  and  four  significant  figures  in  the  final  computed  curve, 
within  the  approximations  inherent  in  the  IBM  program. 

As  a  partial  check  on  the  accuracy  of  the  computing  program,  absolute 
velocities  just  inside  the  free-jet  boundary  at  the  control  points  were 
computed.  If  the  solution  were  exact,  these  velocities  would  all  be  equal 
to  (4e  •  For  H/D  =  1.0,  for  the  five  control  points,  these  velocities 
for  the  converged  solution  varied  from  0.  93  to  0.  99  of  Uc^^with  three  of  the 
five  being  within  3  percent  of  .  For  H/D  ).  Z5,  these  velocities 
varied  from  0.  96  to  1.06  of  ,  with  three  of  the  five  within  3  percent  of 
04,  .  Largest  deviations  from  occurred  for  the  control  points 
closest  to  the  nozzle,  and  hence  closest  to  rt‘gions  where  shape  approxi¬ 
mations  have  been  used.  It  is  believed  that  these  results  could  be 
improved  by  using  more  control  points,  and  by  using  a  finer  breakdown 
for  numerical  integration  purposes.  This  can  be  done  easily,  but  at 
the  expense  of  increased  computer  running  time. 

Another  indication  of  the  accuracy  of  the  analysis  is  provided  by 
Figure  Z8,  in  which  the  absolute  velocity  is  plotted  versus  distance 
from  the  ground  plane  for  H/D  -  1 . 0  at  r/A?  =  1.  7.  The  absolute  velocity 
outside  the  jet  should  be  zero  for  the  correct  jet  boundary  and  correct 
area  vorticity  density  distribution.  For  the  computed  points  shown 
in  Figure  28,  velocities  outside  the  boundary  are  less  than  0.  01  i 
the  calculations  also  show  the  abrupt  change  in  velocity  across  the  jet 
boundary  which  should  be  obtained  theoretically. 

A  detailed  investigation  of  the  various  approximations  in  the  program 
as  they  '•ffect  the  computed  results  was  believed  to  be  unnecessary  for 
the  present  program  in  view  of  the  exi  ellent  agreement  between  computed 
and  experimental  pressure  distributions.  However,  it  is  hoped  th'^t  such 
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an  investigation  will  be  made  because  of  the  potential  application  of  the 
general  technique  employed  in  this  analytical  program  to  other  axially 
symmetric,  three-dimensional,  free-streamline  flows. 

B.  Boundary  Layer  Analysis 

Because  of  the  complexity  of  the  flow  in  the  turning  region  under  the 
impinging  nonuniform  jet,  the  analysis  of  the  flow  along  the  ground  was 
restricted  to  those  regions  where  the  pressure  gradient  at  the  wall 
became  approximately  zero.  Such  a  restriction  required  a  treatment 
which  was  somewhat  different  from  that  presented  for  the  uniform  jet  in 
Reference  1.  The  properties  of  the  wall -jet  portion  of  the  flow  for  both  the 
uniform  jet  and  the  nonuniform  jet  were  computed  using  the  method  of 
analysis  presented  in  this  report.  The  uniform  jet  analysis  predicts 
that  the  wall -jet  velocity  profiles  should  be  given  by  Glauert's  theory  for 
CC  -  1.1  6.  It  was  shown  in  Reference  1  that  the  velocity  profiles 
obtained  with  the  uniform  jet  were  in  excellent  agreement  with  Glauert's 
theory  iord-  1.16,  and  consequently  these  data  will  not  be  reproduced 
herein. 

The  experimental  velocity  profiles  in  the  wall-jet  region  of  the  nonuniform 
jet  with  H/D  =  1.99  are  shown  in  Figure  18  along  with  the  velocity  profiles 
predicted  by  Glauert  for  values  of  cs  =  1.  land  1.2.  These  results  are 
typical  of  the  data  obtained  for  all  values  of  H/D  tested.  The  appropriate 
value  of  a  was  found  to  be  1.  16  in  Section  III  B.  The  accuracy  of  the  pre¬ 
dicted  profile  is  somewhat  less  than  it  was  for  the  uniform  jet.  The 
portion  of  the  data  shown  in  Figure  18  would  appear  to  be  best  represented 
by  a  value  of  ct  slightly  less  than  1.  1  except  in  the  range  >  !• 

The  divergence  from  the  theory  at  large  values  of  was  found  also 

in  the  uniform  jet  work  and  by  other  experiments  and  is  probably  due  to 
a  combination  of  angularity  of  the  flow  with  respi-t  t  to  the  probe  anti  the 
decreasing  accuracy  of  the  measurements  at  the  very  low  velocities 
encountered. 
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The  boundary  layer  region  of  Figure  18  is  plotted  in  Figure  19.  The 
theoretical  profile  (or  a  -^1.  16  has  been  added  to  this  figure.  In  this 
region,  the  experimental  data  are  reasonably  well  represented  by  the 
curve  for  Of  =  1.  16  except  at  values  of  5  1.  33,  where  the  pressure  gradi 
ent  increases  rapidly.  The  curve  for  or  =  1.16  has  been  modified  from 
Glauert's  theory  to  account  for  the  laminar  sublayer.  As  in  Reference  1  , 
the  outer  edge  of  the  laminar  sublayer  is  in  the  region  ^  =^0.  01 S.  The 
theoretical  profile  was  modified  by  matching  the  inner  portion  of  the 
universal  velocity  distribution  for  smooth  pipes  (Figure  20.  4  of  Reference 
5  )  to  Glauert's  profile  at  J/L  =  0.015. 

The  fact  that  slightly  different  values  of  a  appear  to  be  required  in  the 
boundary  layer  and  outer  regions  of  the  velocity  profile  for  agreement 
with  the  experimental  results  may  be  due  to  the  deviations  from  radial 
symmetry  discussed  earlier,  especially  since  the  uniform  jet  tests,  where 
the  radial  symmetry  was  excellent,  provided  extremely  good  agreement 
with  the  theoretical  wall-jet  profiles.  The  small  deviations  from  theory 
in  these  velocity  profiles  constituted  no  serious  problem  in  the  calculation 
of  and  in  Section  111,  since  the  numerical  values  of  the  integrals 
for  nondimensional  mass  flow  and  momentum  used  in  the  analysis  are 
almost  independent  of  cZ  in  the  range  of  Ot  which  is  applicable  to  the 
problem. 

The  calculated  values  of  for  the  uniform  jet  are  compared 

to  the  experimental  data  in  Figure  26(a)  and  26(b)  respectively.  As 
shown,  the  over-all  accuracy  is  quite  satisfactory;  the  difference  between 
theory  and  experiment  never  exi  eeds  10  percent  and  is  generally  closer. 

The  effect  of  varying  H/D  is  to  reduce  the  radius  at  which  the  apparent 
wall-jet  flow  starts.  A  siiinltir  comparison  between  theory  and  experi¬ 
ment  is  made  in  Figure  27(a)  and  27(b)  for  the  nonuniform  jet.  In 
Figure  27(a)  has  bei-n  noncliinensionaliKed  by  and  div  cled  by  2^-  _ 
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The  factor  ]  is  used  only  for  convenience  in  presenting  the  results  on 
logarithmic  graph  paper.  Note  that  the  sample  data  obtained  at  a  second 
value  of  nozzle  mass  flow  (scUd  symbols)  are  essentially  the  same  as  those 
obtained  at  the  mass  flow  used  for  the  bulk  of  the  tests.  Again  the  accuracy 
of  the  approximate  analysis  is  satisfactory  except  near  the  start  of  the 
wall-jet  flow  at  H/D  =  1. 99.  In  this  case  the  assumption  that  the  maxi¬ 
mum  velocity  at  the  start  of  the  wall  jet  can  be  approximated  by  the 
velocity  in  the  free  jet  at  an  equivalent  distance  from  the  nozzle  was  not 
realized,  although  the  use  of  this  assun  ption  provided  reasonable  results 
for  radii  slightly  larger  than  that  corresponding  to  the  start  of  the  wall-jet 
flow.  The  error  in  predicting  for  the  nonuniform  jet  is  somewhat 
larger  than  that  for  the  uniform  jet.  This  is  primarily  due  to  the  higher 
emphasis  placed  on  the  uniform  jet  data  when  estimating  the  mass  flow 
parameter.  Generally,  the  radii  for  which  the  calculated  results  dis¬ 
agreed  most  with  the  experimental  data  coincided  with  the  turning  region 
of  the  jet  as  defined  by  the  ground  static  pressure  distributions  of 
Figure  1  3.  Since  it  was  assumed  in  the  analysis  that  the  jet  was  com¬ 
pletely  turned,  such  results  are  not  unexpected.  Unfortunately,  the 
extent  of  the  turning  region  of  the  non'iniform  jet  at  H/D  -  1, 99  is  larger 
than  the  apparent  radius  for  the  start  of  the  wall  jet,  and  the  calculated 
properties  are  in  error  in  this  region.  The  error  is  on  the  conservative 
side  so  that  the  use  of  the  method  of  analysis  will  at  least  provide  an 
upper  limit  to  the  seriousness  of  the  entrainment  problem  even  in  this 
case. 

It  is  evident  from  the  experimental  results  and  from  the  analysis  that 
the  lack  of  a  method  for  predicting  the  flow  along  the  ground  in  th"* 
turning  region  of  the  jet  is  not  a  serious  limitation.  In  .ill  cases  exce  Pt 
H/D  ^  Z  with  the  nonuniform  jet,  tht“  method  of  analysis  provides  rea¬ 
sonably  accurate  estimates  of  the  properties  jf  the  flow  along  the  ground 
at  the  start  of  th«*  wall-jet  flow.  Sini  e  this  radius  corresponds  to  an 
approximate  maximum  in  velocity  and  a  minimum  in  thickness  of  the  flov/, 
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it  will  in  all  probability,  from  what  is  known  about  particle  entrainment, 
be  the  region  most  critical  to  the  entrainment  problem. 


C.  Comparison  of  Boundary  Layer  Under  Uniform  and  Nonuniform  Jet 

In  Section  III  B. experimental  properties  of  the  flow  along  the  ground 
under  both  jets  were  compared  to  the  results  of  an  approximate 
analysis.  Application  of  this  analysis  requires  a  knowledge  of  the 
conditions  at  the  nozzle  as  a  function  of  H/D  as  well  as  the  free  air 
velocity  decay  curve  o^  each  jet.  However,  the  properties  of  non- 
uniform  jets  at  values  of  H/D  sufficiently  small  that  thf^rt*  is  a  ground 
effect  at  the  nozzle  (H/D.S  2)  are  not  predictable  by  current  theoretical 
techniques.  Hence,  it  would  appear  advisable  to  compare  directly  the 
two  types  of  jet  8  to  discover  if  the  results  from  the  uniform  jet  can  be 
applied  to  the  nonuniform  jet  at  these  low  values  of  H/D.  At  higher  values 
of  H/D,  there  is  no  static  pressure  rise  at  the  nozzle  so  that  these  cases 
require  only  a  knowledge  of  the  free  air  properties  of  the  jet  being  con¬ 
sidered.  The  comparison  will  be  made  on  the  basis  of  constant  thrust 
per  unit  nozzle  area. 


The  above  basis  for  comparison  poses  no  problem  in  the  comparison 
of  for  both  types  of  jets.  The  results  shown  in  Figures  26(b)  and 
27(b)  as  versus  are  independent  of  thrust,  so  they  canbe  com- 

H  r\ 

par ed  directly.  In  the  range  of  H/D  covered  by  the  experiments  and 

r 

for  the  two  sets  of  data  are  essentially  the  same  if  jne  allows  for 

probable  errors  in  the  nonuniform  jet  data  due  to  slight  flow  asymmetry. 
The  results  of  the  two  analyses  are  almost  identical  for  the  cases 
H/  D  =  I  and  2  over  the  complete  range  of  ^  shown.  This  indicates  that, 

as  far  as  the  prediction  of  is  comernid,  the  uniform  j«.’t  analysis 

R 

can  be  applied  directly  to  the  nonuniform  jet  in  ti.o  range  H/D  ^  2.  The 

result  for  the  nonuniforin  jet  will  be  in  error  at  H/D  St  2  for near  the 

R 

start  of  the  wall-jet  flow  for  the  riMsons  discussed  previously  in  the 


31 


TR  64-42 


0 


comparison  of  theory  and  experiment.  At  Jf/D  -  4,  and  in  the  range 
fyH  <  4,  both  the  experimental  and  analytical  results  indicate  that 
the  two  jets  arc  not  comparable,  the  flow  under  the  uniform  jet  being 
thinner  than  that  under  the  nonuniform  jet.  Hence,  for  H/D  >  2,  the 
decay  in  maximum  velocity  in  the  nonuniform  jet  must  be  accounted 
for  if  reasonably  accurate  predictions  are  to  be  obtained. 


In  order  to  compare  the  maximum  velocity  near  the  ground  generated 
by  the  two  jets,  it  is  necessary  first  to  convert  the  results  presented 
in  Figures  26(a)  and  26(b)  to  ,  where  r  is  the  thrust  and  A  is 

nozzle  area  (n^*)  .  For  the  uniform  jet,  this  is  done  by  expressing 
the  thrust  in  terms  of  .  It  can  be  shown  that 


and  hence 


1'-(t 

^  C  ^  *^.1  ff  / 


)‘ 


The  variation  of  —  with  H/D  is  shown  in  Figure  29  ,  for  a  uniform 
^  /  /  / 

The  limit  of  — ^  as  H/D  approaches  zero  was  calculated  from 


jet. 


considerations  of  two-dimensional  inviscid  flow  through  a  sluice  gate. 
Since  reached  a  maximum  value  of  approximately  1  at  the  start 

of  the  wall  jet  for  all  H/D<  4,  the  maximum  veloc  ity  reached  at  the 
ground  is  simply  u._  / - - — .. 

for  the  uniform  jet.  The  value  of  shown  in  Figure  30 

for  both  the  uniform  jet  and  the  nonuniform  jet.  The  points  for  the 


latter  were  taken  directly  from  Figure  27(a)  because*  the  experimental 

data  for  the  nonuniform  jet  indicated  that,  in  the  range  H/D—  ]. 

*  Figure  29  may  also  be  used  tc  obtain  an  c*stimate  of  the  thrust 
variation  with  constant  power  input  for  an  axially  symmetric, 
uniform,  impinging  jet.  Using  Fquations  (14)  and  (15)  .ind  the  recjuire 
ment  that  power  remain  constant,  one  can  show  1  .  t^/Q.  1^1 

1\  ij  )  1'*(0J  J 

for  constant  power.  This  equation  plus  Figure  29  may  be  used  to 
calculate  the  variation  of  as  a  function  of  H/D. 
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TTp/f^U^^-  ^  and  hence  •  The  point*  shown  are 

the  calculated  velocity  at  the  start  of  the  wall  jet  (break  in  the  curves  of 
Figure  ZTl.i)  ), 

It  can  be  seen  that  the  uniform  jet  analysis  provides  results  for 

identical  to  those  obtained  with  the  nonuniform  jet  analysis  in 
the  range  H/Di!  2.  In  view  of  this  result  and  the  corresponding  result;* 
for  t  it  can  be  assumed  with  reasonable  certainty  that  the 

uniform  jet  analysis  can  be  applied  to  calculate  the  properties  of  the 
flow  along  the  ground  under  nonuniform  jets  for  values  of  H/D  <  2. 

The  results  of  the  analysis  will  of  course  be  valid  only  for  radii  suffi¬ 
ciently  large  that  the  flow  has  been  completely  turned.  At  values  of 
H/D  £2,  the  velocity  decay  characteristic  of  the  nonuniform  jet  ex¬ 
hausting  into  free  air  must  be  included  in  the  analysis. 

A  general  comparison  of  both  the  analytical  and  the  experimental  results 
for  the  two  jets  provides  the  following  features.  At  values  of  H/D  S  !• 
the  flows  along  the  ground  under  the  two  jets  are  essentially  the  same. 
The  limit  H/D  ^  1  rather  than  |  if>  stipulated, even  though  quantitative 
measurements  of  the  flow  along  the  ground  were  not  obtained  with 
the  nomniform  jet  at  H/D  =  1.  Flow  visualization  pictures  and  the 
ground  static  pressure  distribution  obtained  at  H/D  =  1  indicated  that 
the  effective  turning  radius  of  the  jet  was  sufficiently  small  as  to  caus  ; 
no  such  problems  as  those  observed  at  H/D  =  2.  At  values  of  H/D>  1, 
the  flow  along  the  ground  under  the  nonuniform  jet  had  a  lower  maximum 
velocity  and  larger  thickness  than  that  under  the  uniform  jet.  These 
differences  were  not  large,  but  they  do  provide  the  unexpected  result 
that  the  nonuniform  jet  is  slightly  less  prone  to  cause  particle  entrain¬ 
ment  at  large  values  of  H/D,  and  ii  no  worse  than  the  uniform  jet 
at  small  values  of  H/D. 


i 
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APPENDIX  I 


GENERAL  COMMENTS 
ON  THE  FLOW  VISUALIZATION  PICTURES 


Splitter  Plate  Pictures 

The  streamline  patterns  obtained  on  the  splitter  plates  must  not  be 
taken  as  an  exact  picture  of  the  flow  which  wou^d  occur  in  the  absence 
of  such  a  plate.  Since  a  surface  has  been  inserted  in  the  flow,  there 
will  be  a  boundary  layer  developed  on  this  surface.  In  general,  the 
combined  thickness  of  the  plate  and  the  displacement  thicknesses  of 
the  boundary  layers  developed  on  each  side  of  the  plate  will  be  a 
negligible  fraction  of  the  jet  diameter. so  that  this  thickness  will  have 
little  effect  on  the  flow  outside  the  boundary  layer.  However,  in 
regions  where  there  is  a  pressure  gradient  in  the  flow  (close  tc  the 
ground  in  this  cast*),  the  off'  !  vhe  pressure  gradient  on  the  flow  in 

e  boundary  'aver  will  be  t  celerale  or  decelerate  mostly  that 

pc’^tion  of  Iht*  air  with  th<‘  iec  it  kinetic  energy.  This  is,  of  course, 
the  region  adjacent  to  the  surface,  the  one  instrumental  in  producing 
the  streamline  patterns 

The  net  effect  of  the  pressure  gradients  in  the  flow  near  the  ground 
will  be  twofold.  First,  near  the  jet  centerline  the  adverse  pressure 
gradient  occurring  as  the  ground  is  approached  will  separate-  the 
boundary  layer.  Such  behavior  was  obsc-rved  in  the  uniform  jet  tests 
and  is  the  reason  why  t!ie  trailing  edge  of  the  splitter  plate  was  kept 
at  least  6  inches  from  the  ground  b.  ar<i.  In  this  case  the  flow  in  the 
core  of  the  jet  ’.vas  essentially  free-  of  turbulence  and  the  velot  ity  at 
whi'  h  the  tests  were  conducted  wa.s  sufficiently  low  so  as  to  sustain 
a  lamin.'.r  boundary  layer  near  the  jc-t  axis  .ind  hence  to  provi  de  a  high 
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susceptibility  to  separation.  With  the  nonuniform  jet,  the  turbulence 
in  the  flow  probably  caused  a  forced  transition  to  a  turbulent  boundary 
layer  near  the  plate  leading  edge.  The  higher  resistance  of  a  turbulent 
boundary  layer  to  separation  and  the  less  adverse  pressure  gradient;! 
observed  with  this  jet  explain  the  absence  of  separation  in  the  flow 
pictures  which  were  made  with  the  trailing  edge  of  the  splitter  plate 
only  1  inch  from  the  ground  board. 

The  second  effect  is  that  of  the  pressure  gradient  in  the  radial  direction 
in  the  absence  of  separated  flow.  In  this  case,  the  flow  in  the  boundary 
layer  very  near  the  wall  will  be  turned  more  than  that  in  the  free  stream, 
again  because  of  the  difference  in  kinetic  energy.  This  will  show  up  in 
the  flow  pictures  as  an  exaggerated  turning  of  the  flow  near  the  ground 
and  possibly  an  increased  apparent  thickness  of  the  radial  outflow  along 
the  ground.  The  portions  of  the  flow  patterns  ir  I’le  center  of  the 
impingement  region  under  the  nonuniform  jet  shoulc'  be  least  affected 
by  the  interaction  of  boundary  layer  and  radial  pressure  gradient,  as 
this  region  appeared  to  be  at  roughly  constant  pressure. 

The  above  comments  apply  only  to  the  patterns  obtained  on  the  splitter 
plates.  The  conclusions  to  be  drawn  are  that  the  patterns  must  be 
regarded  as  providing  only  an  approximate  picture  of  the  flow  in  any 
region  where  the  streamline  patterns  show  a  large  curvature,  with  the 
possible  exception  ol  flow  in  the  center  of  the  impingement  area  under 
the  nonuniform  jet. 

Ground  Flow  Pictures 

The  streamline  patterns  shown  reflect  the  flow  in  the  boundary  layer 
along  the  ground.  Since  the  surface  used  to  obtain  the  patterns  is  in  the 
flow  at  all  times  and  is  not  an  additional  surface  as  with  the  splitter 
plate,  the  resulting  streamline  pattern  is  reasonably  representative  of 
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the  actual  flow,  except  for  a  gravity  effect.  The  board  was  mounted 
vertically  in  these  tests  ;  hence,  in  regions  where  the  flow  velocity 
was  low,  gravity  caused  the  lampblack  and  kerosene  mixture  to  flow 
downward.  This  is  especially  apparent  in  the  center  of  the  impingenr. ent 
region  where,  in  some  cases,  there  is  a  pronounced  downward  curvature 
to  all  of  the  streamlines. 
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APFKNDiX  II 


DETAILS  OF  UNIFORM  INVISCID  IMPINGING  JET  ANALYSIS 


Derivations  of  Equations  for  Inviscid  Impinging 
Jet  FIonv  Basic  Formulation 


Consider  a  Cartesian  coordinate  system  with  the  z-axis  directed 
along  the  axis  of  symmetry  of  the  flow,  positive  vertically  up 
from  the  ground  plane,  and  the  x-y  plane  (  z  -  0)  coincident  with 
the  ground  plane.  Consider  an  eiement  of  volume 
with  vorticity  vector  i').  For  vorticity  distributed  through 

out  a  given  /olume  /  ,  one  can  write  the  vector  velocity,  ^  , 

at  any  point  (  x,  y,  z)  in  the  fluid  as  (RefiTem  e  6) 


FP 


(ii 


where 


d  -  (z  -  x')  0  *  ( i-  ( i  - 


AAA 

and  t.y.k  are  unit  vectors  in  the  x,  y,  and  z  dirt'itions,  respect 
ively.  In  the-  case  of  .ixially  svniinetric  flow, 

j  ~  C*,  ^  Cy  J 

The  trans  tor  in.iticjn  to  .  ylindrical  coordinates  {  f,  i  ,  S  )  where 

X  •  r  cos  ^  dLr\  0 

a  nd 

Cjc  *  Co  ^  •  Cu  ~  •  Co  ^ 


gives 


2  rr'cos&'y^* 

where  t^,j0  and  are  unit  vectors  in  the  ^  Q  and  2 

directions,  respectively.  Because  of  axial  symmetry,  only  the 
6*0  plane  need  be  considered. 


Substituting  in  Equation  (II- 1)  and  equating  vector  components. 


we  have 


um-—f  Hr- f')  COS  9'dr(r',§\  G') 

-  Zrr'cosO']^* 

J_  f  i')(rcosO'-r')dv(r'i\  9') 


(II-2) 


where  u, ,  v*  are  the  velocity  compone  ts  parallel  to  the  r-axis 
and  the  z-axis,  respectively,  and  is  the  velocity  component 
perpendicular  to  the  meridian  plane.  The  integration  is  taken  over 
the  volume  which  includes  all  the  distributed  vorticity  of  intensity 
<o  (r,l). 


Consider  now  the  vorticity  distributed  in  a  thin  layer  on  the  edge  of 
the  jet.  Consider  a  given  element  oi  volume,  ,  of  this  layer 

(Figure  31)  whose  projection  on  the  tneiudian  plane  is  €CL5  ,  where 
€  IS  the  infinitesimal  thickness  ol  the  l^yer  and  cis  if  the  length 
of  the  element  along  the  curve  r*  ^q(i)  which  defines  the  edge  of 
the  jet.  From  Stoke's  theorem  (Reterence  6,  page  46)  for  a  circuit 
C  in  the  meridian  plane  enc  losing  the  projection  of  dt  t  we  have 

f  n-CdA  -  f  q  -db 


7 

Jc 


where  fi  is  the  unit  outward  normal  vector  to  the  surface  enclosed 
by  C,  s'  is  the  area  enclosed  by  C’  and  ^  is  the  velocity  vector 
around  C  .  The  integral  on  the  right-hcnd  side  is  taken  in  the 
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I 


clockwise  direction.  Then 

I  ri'^dLA  =  (^^edL^ 

Js' 

hence,  ^ 

where  and  arc  as  shown  in  Figure  31.  If  we  let  € -^  Ot 
such  that 

3  constant 

then 


Also, 


^^er^dsdff' 


dide'  ■ 


Thus,  Equation  (11*2)  can  be  written 


f 

L  di  J 

\z-e)co5Q'di  to* 

*^[r0(i')y~  Zrrf^(f)  cos9'] 

17i 

(II-3) 

dl'd9' 

H 

dtod') 

[rcoiS'-r^Cz')]^ 

*  -  2rrf,(t')cos9'[ 

pTiC- 

where  the  integration  is  carried  out  over  the  entire  vortex  sheet  $ 


The  Boundary  Conditions  for  Infinite  Impinging  Jet 

The  boundary  conditions  on  the  inviscid  impinging  jet  flow  are  as  follows; 
Along  the  surface, z  =  0,  I/"  -  0.  Along  the  jet  boundary  which  is  a 
free  streamline,  the  velocity  is  equal  to  a  constant:  therefore, 

/  is  also  constant,  as  just  outside  the  jet  the  velocity  is  zero. 
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To  satiffy  the  boundary  condition  that  ^  =  0  for  z  =  0,  the  plane 
z  =  0  it  made  a  plane  of  symmetry  for  the  impinging  jet  and  its 
image.  Hence,  the  integration  over  j  ,  Equation  (11-3).  must  be 
carried  out  over  two  vortex  sheets,  that  defining  the  impinging 
jet  and  that  for  the  image  jet. 

The  Boundary  Conditons  for  Jet  Nozzle  at  Finite  Height 
Above  the  Ground 

The  jet  is  assumed  to  issue  from  a  straight  tube  of  constant 
diameter  (Figure  32)  at  finite  height,  H  ,  above  the  ground.  A 
free-stream  surface  issues  from  the  edge  of  the  jet  exit,  B  in 
Figure  32,  such  that  the  free  streamline  velocity  vector  at  B  is 
parallel  to  the  side  of  the  tube.  The  flow  inside  the  tube  along  thd 
side  of  the  tube  must  be  tangential  to  the  side;  but,  of  course,  the 
streamline  along  this  side  is  no  longer  a  free  streamline,  as  a 
pressure  difference  can  exist  and,  hence,  the  magnitude  of  the 
edge  velocity  may  vary.  The  boundary  condition  along  the  e.?<(e  of 
the  tube  is,  then,  that  the  normal  component  of  velocity,  ca  ,  is 
zero.  In  addition,  at  B  in  Figure  32,  t/'  is  continuous.  Again, 
along  the  free  streamline  BC,  the  velocity  and,  hence,  ^  are  con* 
stant. 

Final  Form  for  Velocity  Components 

Figure  32  illustrates  the  flow  model  as  thus  far  developed.  The 
straight  tube  of  constant  diameter,  D  ,  centered  on  the  z-axis 
extends  from  z  =  H  to  z  -  o»  .  For  0<l<  H  ,  the  curve  t  • 
defines  the  edge  of  the  jet.  To  insure  that  v*  is  tangential  at 
B  i.'>  Figure  32,  we  require  that  *0  at  z  =  H;  that  is,  there  is 
no  discontinuity  in  slope  between  the  side  of  the  no7zle  and  the  free* 
stream  surface.  The  sides  of  the  tube  and  the  jet  are  represented 
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by  a  vortex  sheet  of  vorticity  ^(Tq,z)  per  unit  area  whose  r-  and  z- 
components  are  zero.  The  image  of  this  vortex  system  for 
assures  that  the  boundary  condition  (no  normal  flow)  along  OD 
is  satisfied. 


The  integrations  with  respect  to  ^^can  be  done  in  closed  form  in  each 
of  the  integrals  of  Equation  (11-3).  If  this  is  done,  and  we  let  , 

and  there  results 


^  .1  1.., 

•  ■*  IFfU  v:— v.^  ^ _ . - ,  ff*, )-*■«,  d  i 

Jo  /Vs 


df 


[TpiriX  . 

j  ,  f=  =  .=/xU 


Ji 


-/ 


E(k4)-K(k^) 


V(F,V 


^  iJf  i(F  -  (r*  f?  [  if-  F')  *-(f  -  f)  J 


-t 


HT’rTTMl^ 


(i-r)^*rHro(r)]‘  . 


d/' 


H- 


redro(f')] 

df  J 


£(k^)-<(A^)ldg'}  (11-4) 
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where  K(k)  i*  the  complete  elliptic  integral  of  the  first  kind  and 
£(k)  is  the  complete  elliptic  integral  of  the  second  kind,  with  modu¬ 
lus  k  ;  also. 

n  ^  tk  I* 

'(TTFpTTTTTp 

/L  »I  (k  t» 

These  integrals  are  valid  everywhere  except  on  the  boundaries,  at 
(l,f)  for  and  (F^(i\l)ioT  / 5  *  where  the  integrands  are 

singular.  However,  even  on  the  boundaries,  it  can  be  shown  that  the 
Cauchy  principal  values  of  these  integrals  exist. 

Equations  (11-4)  enable  the  velocity  at  any  point  in  the  flow  to  be  com¬ 
puted,  once  the  form  of  fgd) ,  the  curve  defining  the  boundary  of  the 
jet,  and  ?(f)  ,  the  vorticity  area  density  of  the  vortex  sheet  on  the 
boundary  of  the  jet  and  the  tube,  are  known. 

Equations  (11-4)  with  the  as'ociated  boundary  conditions 

Um  u,(r.I)-0,  t±HlR»  f<1 

F-^f 


a(F,r)  dF^(%) 

“  dl  * 


F<ro 


give  rise  to  simultaneous  nonlinear  .  utegral  equations  for  ^(fl)and 

,  which,  as  they  stand,  do  not  admit  of  solution  in  closed  form. 
However,  by  means  of  a  high-speed,  large-capacity  digital  computer, 
a  solution  might  be  possible  if  a  convergent  iterative  process  can  be 
derived  which  will  permit  the  determination  of  the  jet  boundary  . 
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B.  Detailw  of  IBM  704  Computer  Program 


The  general  computational  procedure  adopted  it  at  lollowt: 

1.  An  initial  shape  of  the  free-ttreamline  boundary  it  attumed. 

2.  The  boundary  condition  that  the  flow  at  the  wall  of  the  tube 
is  parallel  to  the  tube  it  uted  to  determine  the  varying  area 
vorticity  density  of  the  vortex  theet  representing  the  tube 
and  the  constant  area  vorticity  den^^ity  of  the  free>ttream 
surface  vortex  sheet. 

3.  The  flow  normal  to  the  free-s^ream  surface  is  then  evaiuateoi 
if  there  is  a  velocity  normal  tc  the  boundary  assumed  for 

the  tree-stream  surface,  that  assumed  boundary  is  incorrect 
and  must  be  adjusted. 

4.  The  iteration  proceeds  for  ihe  adjusted  boundary  (  by  going 
back  to  step  1). 

The  boundary  condition  at  the  wall  of  the  tube,  step  2,  it  satis¬ 
fied  at  a  discrete  number  of  points,  or  values  of  J  ,  along  the 
tube.  The  nondimens’onal  area  vorticity  oensity,  id)  in  Equation 
(11-4),  is  assumed  to  be  a  constant  within  the  increment  cen¬ 
tered  on  each  of  these  points  (but,  of  course  varying  from  point 
to  point).  The  firs*  two  integrals  in  the  expre  ssion  for  V,  Equa¬ 
tions  (11-4),  can  then  be  expressed  as  a  sum.  Setting  V-  0  at  each 
of  the  specified  values  of  \  ,  there  results  a  series  of  linear,  simul¬ 
taneous  equations  in  the  unknown  id) .  For  f large  enough,  id) 
approaches  a  constant  limiting  value,  and  this  fa^'t  is  utilized  to 
limit  the  number  of  simultaneous  equations  which  must  be  solved. 

It  happens  that  if  fd)  is  a  constant,  the  first  two  integrals  in  the 
equation  for  tA  in  Equations  (11  4)  can  be  integrated  -explicitly,  which 
eases  the  numerical  integration  process  on  the  computer  and  also 
lends  considerably  to  the  accuracy  with  which  IT  c.in  be  computed. 
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f 


f 


V 


In  practic  e,  the  value  of  /  at  which  approaches  a  constant 

is  initially  guessed  at,  and  then  adjusted  after  initial  runs  have  been 
made.  The  computed  distributions  of  do  not  appear  to  vary 

significantly  with  variations  in  the  assumed  jet  boundary  for  a  given 
H/D,  at  least  for  the  two  cases  computed  in  the  present  work. 


The  process  of  determining  the  revised  jet  boundary  in  step  3  is 
based  on  the  boundary  condition  along  the  jet  boundary  free  stream  • 
lines,  that  the  velocity  normal  to  the  boundary  is  zero.  For  example, 
if  the  computed  velocity  normal  to  the  jet  boundary  is  directed  out¬ 
ward,  this  intuitively  would  suggest  that  the  assumed  boundary  should 
be  displaced  outward.  Let  the  curve  for  the  assumed  free  jet  boundary 
BC  (Figure  32)  for  the  i^^  repetition  of  steps  1  through  3  be 

(^o)i  • 

Once  the  distribution  of  the  vorticity  in  the  vortex  sheet  is  known 
from  step  2,  above,  it  is  possible  to  compute  the  velocity  vector 
anywhere  in  the  flow,  based  on  ,  the  assumed  free-stream 

surface  boundary  using  Equations  (iI-4).  The  velocity  components 
U(F^,B)  ,  at  the  boundary  are  computed;  then 

the  slope  of  the  velocity  vector  along  the 
assumed  boundary.  If  this  slope  is  equal  to  d(F^)^/cLf  all  along 
the  assumed  boundary,  is  the  correct  (  urve,  as  all  boundary 

conditions  are  then  satisfied.  If  not,  we  let 


d 


h(f)i 


3 


(II-5) 


where  /C  is  a  constant,  called  the  "gain"  constant.  It  should  be 
noted  that  no  mathematical  rigor  is  claimed  for  this  procedure;  it 
is  primarily  a  s ystemization  of  the  method  which  a  computer  might 
use  to  adjust  boundaries  in  a  successive  approximation  technique. 
Convergence  can,  thus  far,  be  determined  only  by  a  trial  and  error 
process.  Several  iterations  are  made  with  a  selected  value  of  the 
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constant.  If  the  successive  boundaries  appear  to  be  diverging,  then 
successively  smaller  values  of  ^  are  used  until  convergence  is  ob¬ 
tained. 

Although  the  integrations  indicated  in  Equations  (11-4)  are  all  taken 
with  respect  to  the  coordinate  i  ,  in  actuality  the  integrals  with 
limits  of  integration  from  0  to  H/R  were  broken  up  into  several  integrals; 
see  Figure  33.  In  region  I,  the  integration  is  with  respect  to  2 
between  the  limits  fo  and  H/R.  For  regions  11,  III,  and  IV,  the  in¬ 
tegral  is  transformed  so  that  integration  is  with  respect  to  f'  ,  with 
limits  as  shown  in  Figure  33.  In  each  of  the  noted  regions,  breakdowns 
for  purposes  of  numerical  integration  can  be  varied  individually. 

Certain  assumptions  and  approximations  are  involved  in  the  actual 
computer  program.  Equation  (11-5)  is  applied  at  each  iteration  at 
specified  values  of  ^  ,  called  "control  points,  "  in  regions  II  and  III 

in  Figure  33,  These  control  points  may  or  may  not  coincide  with  values 
of  A  used  in  numerical  integrations.  Once  the  jet  boundary  at  the 
control  points  has  been  adjusted  in  accordance  with  Equation  (11-5), 
increments  are  added  to  the  jet  boundary  curve;  these  increments 
vary  linearly  between  the  control  points  and  match  the  increments  at 
the  control  points  computed  from  Equation  (11-5).  The  value  of  is 
kept  as  close  to  unity  as  possible,  consistent  with  accuracy  of  integra¬ 
tion.  Variation  of  (i)  between  F  -  ^nd  f  •  H//f  is  based 

on  the  variation  of  at  the  control  point  closest  to  ,  although 

the  shape  is  determined  by  the  initial  assumed  curve.  For  ?  >f^  , 

the  shape  of  is  assumed  to  vary  as  C/r  ,  with  C  a  con¬ 
stant,  determined  from  Equation  (14).  Integration  is  cut  off  at  ; 

the  value  of  is  made  as  large  as  necessary  so  as  not  to  affect  the 
accuracy  of  velocities  computed  at  the  control  points.  The  computations 
for  which  results  are  presented  in  this  report  were  made  with  five  control 
points  ,  For  H/D=  i.  0  ,  r^-X.Z  aixl  =  4.  4;  for  H/D  -  0,  Z5,  r^=  1.05 and 

=  4.  5.  In  both  cases,  Fy  =  50. 
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C.  Outline  of  a  Po«»iblc  Improved  Iteration  Techniqute 


uCr^  ,i)  -  yi)  *0 


Consider  the  flow  tangency  boundary  condition  on  the  free  streamline 

(U-6) 

For  an  assumed  free>streamline  curve  (.^9)0,  other  than  that  for  the 
exact  solution,  there  would  result 

<U-7) 


‘£(1). 

For  purposes  of  numerical  integration  on  a  digital  computer,  the 
integrals  constituting  u-  and  I/*  are  evaluated  by  computing  the 
integrand  at  specific  values  of  /  ,  say,  ,  and  summing  th^ 

results  by  an  appropriate  rule  such  as  Simpson's  rule.  On  this  basis. 
Equation  (11-7)  can  be  stated  (A^  i"  •  if^) 

**  . 

£u.J(F,)^i  i  (r,)« ,  t„  )  (n-8) 

M  m  ^  * 


i-0 


^(ro) 


o  /an 


df 


(11-9) 


*•0 


Now,  let 

f  ury)i 

where  (^o)i  is  the  exact  value  of  and 

is  the  incremental  error  in  at  §  .  If  Equation  (II-9) 

is  substituted  in  Equation  (11-8)  and  the  terms  ,  ^1/'^*  are  expanded 

to  first  order  in  ,  there  results 

<-*d  ^ 


4^0 


«  £n(ln)  , 


for  n*  0  to  M 

(U-IO) 
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where  equ.ition  (ll-l)  has  been  applied;  Six  I  >  3re  the 

eoeffic ients  of  the  first-or  ler  terms  in  of  the  expansions  of 

Equation  (II -10)  for  r7»^to  n  *  JD  results  in  M  +  1 
linear,  simultaneous  algeb  raic  equations  in  the  M  1  unknowns 

;  i  i  '0  to  M.  However,  the  (’^j^-are  also  unknown;  hence, 
ihe  substitution 

is  made  in  Equation  (II- 10), 

This  is  the  basis  of  the  proposed  iteration  scheme.  The  initial  assumed 
shape  of  the  free-streamline  boundary  [(^o  )a,i  ] ^  ts  used  in  Equation 
(11-7)  to  compute  ,  and  the  M  +  1  equations  resulting  from 

Equation  (11-10)  (with  Equation  (II- 11)  included)  are  solved  for  the 
M  +  1  f(Ar^  .  The  process  is  repeated  v/ith 

The  iteration  proceeds  until  all  the  ];  are  as  small  as  desired. 

Whether  convergence  is  obtained  at  all  undoubtedly  depends  on  how 
close  the  assumed  boundary  is  to  the  exa^t  boundary.  However,  there 
seems  to  be  no  quest  on  that  if  convergence  were  obtained  by  the 
iteration  method  actually  used  during  the  present  calculations  (empiric¬ 
ally  determined  "gain  constant,"  Equation  (II  -  5)  ), convergence  would  be 
obtained  much  more  rapidly  by  the  method  proposed  here.  It  is  noted 
that  if  all  fAe, except  )„  had  been  assumed  to  be  zero,  then  an 

analytic  expression  for  the  empirical  "gain  constant"  (Equation  (II-5)  ) 

would  result.  Hence,  such  a  value  of  at  each  would  undoubt¬ 

edly  be  an  optimum  value  for  most  rapid  convergence  using  Equation 
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(II-5).  However,  the  proposed  method  based  on  Lquations  (U-IO) 
and  (II-ll'  should  result  in  even  faster  convergence.  Unfortunately, 
this  method  was  devised  too  la’^e  in  the  present  research  to  be  in¬ 
corporated  in  the  IBM  704  proprnm.  It  is  hof*"d  that  in  the  near 
future,  there  will  be  opportunity  to  employ  the  method  in  the 
solution  of  a  different  problem. 
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3268 


Figure  I  SHEAR  SCREEN  ASSEMBLY 
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Figure  2  FLOW  VISUALIZATION  APPARATUS 
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Figure  4  FLOW  VISUALIZATION  PICTURES  OF  NORMALLY  IMPINGING  NONUNIFORM  JET,  H/D  =  1.99 

(b)  FLOW  ALONG  GROUND 
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Figure  5  FLOW  VISUALIZATION  PICTURES  OF  TILTED  IMPINGING  NONUNIFORN  JET,  H/O  =  1.99 

(a)  FLOW  IN  PLANE  OF  SYMMETRY 
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Figure  5  FLOW  VISUALIZATION  PICTURES  OF  TILTED  IMPINGING  NONUNIFORH  JET,  H/D  =  1.99 

(b)  FLOW  ALONG  GROUND 
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Figure  7  VaOCITY  PROFILES  IN  FREE  NONUNIFORM  JET  AT  NOZZLE  EXIT 
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Figure  10  VARIATION  OF  HASS  FLOW  IN  UNIFORM  AND  NONUNIFORM 
FREE  JETS  WITK  DISTANCE  FROM  NOZZLE 


61 


TR  6-4-42 


Flgurtll  MEAN  VELOCITY  PROFILES  IN  NOZZLE  EXIT  PLANE  FOR  VARIOUS  H/D. 
NONUNIFORM  JET 
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SYNML 


Fi«ur«  13  GKOUNO  STATIC  PtESSWE  DISTtltUTIONS,  NOtMAUY  IMPIlieiM  HOMWIFOtM  JH 

(b)  H/0  =  0.99 


Figure  13  6tOUIID  STATIC  PRESSUtE  DISTRIBUTIONS,  NORMALLY  IMPINGING  NONUNIFORM  JET 


Figure  IS  STATIC  PtCSSUIE  DISTRIBUTIONS  ON  JET  CENTERLINE,  NONUNIFORN  JET 


igur*  16  AISOLUTE  VaOCITY  DISTRIBUTIONS  IN  TURNING  REGION  OF  IMflNGING  NONUNIFORM 
JET  AT  SEVERAL  CONSTANT  DISTANCES  FROM  GROUND,  H/D  =  1.99 
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Figurt  17  VELOCITY  FKOFILES  NEAI  MOUND,  NOtNALlY  INFINOINQ  NONUNIFOIM  JET 

(a)  N/D  X  0  00 
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Figurt  17  VELOCITY  PROFILES  *(EAR  GROUND,  NORMALLY  IMPINQINO  NONUNIFORN  JET 

(b)  H/0  s  1.99 
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Figurt  17  VELOCITY  PROFILES  NEAR  OROUND,  NORMALLY  INPINOINO  NONUNIFORM  JET 

(c)  H/D  =  N.07 
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Figur*  21  JET-CENTERLINE  STATIC  PRESSURE  DISTRIIUTIONS,  UNIFORM  JET 
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Figure  25  NO  NO  I  HENS  I  ON  AL  RADIAL  MASS  FLOM  vt.  r/R  FOR  lOTH 
UNIFORM  AND  NONUHIFORH  JETS 
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Figure  26  VARIATION  OF  FLOM  PROPERTIES  NEAR  GROUND  WITH  RADIAL  DISTANCE 
FROM  STAGNATION  POINT,  NORMALLY  IMPINGING  UNIFORM  JET 
(•)  MAXIMUM  VELOCITY  NEAR  GROUND 


VARIATION  OF  FLOW  PROFERTIES  NEAR  OROUNO  WITN  RADIAL  DISTANCE 
FROM  STAONATION  POINT  NORMALLY  IMPIN6IN6  UNIFORM  JET 
(b)  FLOW  THICKNESS 
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Figure  27  VARIATION  OF  FLOW  PROPERTIES  NEAR  GROUND  WITH  RADIAL  DISTANCE 
FROM  STAGNATION  POINT,  NORMALLY  IMPINGING  NONUNIFORM  JET 
(e)  MAXIMUM  VELOCITY  NEAR  GROUND 
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Figure  27  VARIATION  OF  FLOW  PROPERTIES  NEAR  GROUND  WITH  RADIAL  DISTANCE 
FROM  STAGNATION  POINT,  NORMALLY  IMPINGING  NONUNIFORM  JET 

(b)  FLOW  THICKNESS 
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VARIATION  OF  U/U^  WITH  H/O  FOR  A  UNIFORM  JET 
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Figure  32  FLOW  HO  Da  USED  IN  INVISCIO  ANALYSIS  OF  CIRCULAR 
IMPINGING  UNIFORM  JET 
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U.  S.  Army  Combat  Developments  Command 

Aviation  Agency  1 

U.  S.  Army  War  College  1 

U.  S.  Army  Command  and  General  Staff  College  1 

U,  S.  Army  Aviation  School  1 

U.  S.  Army  Infantry  Center  Z 

U.  S.  Army  Aviation  Test  Board  1 

Air  Force  Systems  Command,  Wright-Patter son  AFB  1 

Air  Force  Flight  Test  Center,  Edwards  AFB  1 

Bureau  of  Naval  Weapons  4 

U.  S.  Naval  Postgraduate  School  1 

Naval  Air  Test  Center  1 

David  Taylor  Model  Basin  1 

Ames  Research  Center,  NASA  1 

NASA-LRC,  Langley  Station  1 

Lewis  Research  Center,  NASA  1 

Manned  Spacecraft  Center,  NASA  1 

NASA  Representative,  Scientific  and  Technical 

Information  Facility  Z 

Research  Analysis  Corporation  1 

National  Aviation  Facilities  Experimental  Center  1 

U.  S.  Army  Standardization  Group,  Canada  1 

Canadian  Liaison  Officer, 

U.  S.  Army  Transportation  School  1 

British  Army  Staff,  British  Embassy  1 

U,  S.  Army  Standardization  Group,  U.  K.  1 

Defense  Documentation  Center  10 
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